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[l]  Interplanetary  coronal  mass  ejections  (ICMEs)  are  the  primary  cause  of  severe  space  weather  at 
Earth  because  they  drive  shocks  and  trigger  geomagnetic  storms  that  can  damage  spacecraft  and 
ground-based  systems.  The  Solar  Mass  Ejection  Imager  (SMEI)  is  a  U.  S.  Air  Force  experiment  with  the 
ability  to  track  ICMEs  in  white  light  from  near  the  Sun  to  Earth  and  beyond,  thus  providing  an  extended 
observational  range  for  forecasting  storms.  We  summarize  several  studies  of  SMETs  detection  and 
tracking  capability,  especially  of  the  ICMEs  associated  with  the  intense  (peak  Dst  <  -100  nT)  geomagnetic 
storms  that  were  the  focus  of  the  NASA  Living  With  a  Star  Geostorm  Coordinated  Data  Analysis 
Workshop.  We  describe  the  SMEI  observations  and  analyses  for  the  18  intense  storms  observed  from  May 
2003-2007  with  adequate  SMEI  coverage  and  identified  solar  and  interplanetary  source  regions.  SMEI 
observed  the  associated  ICMEs  for  89%  of  these  intense  storms.  For  each  event  we  extracted  the  time 
differences  between  these  sets  of  times  at  1  AU  for  shock  arrival  time,  predicted  ICME  arrival  time,  onset  of 
high-altitude  aurora  observed  by  SMEI,  and  storm  onset.  The  mean  intervals  between  successive  pairs 
of  these  data  were  found  to  each  be  hours.  On  average,  SMEI  first  detected  the  geoeffective  ICME 
about  1  day  in  advance,  yielding  a  prediction  lead  time  of  ~18  hours.  Finally,  the  RMS  values  for  the 
ICME-shock  and  storm-ICME  time  differences  were  determined,  and  provide  at  least  a  1-hour 
improvement  compared  to  similar  observational  and  model-dependent  studies. 
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1.  Introduction 

[2]  Coronal  mass  ejections  (CMEs),  which  are  termed 
interplanetary  coronal  mass  ejections  (ICMEs)  when  they 
traverse  interplanetary  (IP)  space,  are  the  primary  source 
of  severe  space  weather  at  Earth.  They  cause  geomagnetic 
storms  that  can  damage  both  military  and  civilian,  space 
and  ground  assets.  These  hazardous  storms  are  difficult  to 
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forecast  and  there  are  many  false  alarms  [e.g.,  McKcnna- 
Lawlor  et  al.f  2006].  The  Solar  Mass  Ejection  Imager  (SMEI) 
experiment  [Eyles  et  al,  2003;  Jackson  et  al ,  2004]  on  the 
Coriolis  spacecraft  has  been  obtaining  white  light  images 
of  nearly  the  full  sky  every  102  min  for  the  last  5  years. 
SMEI  views  the  sky  above  Earth  using  sunlight-rejecting 
baffles  and  charge-coupled  device  (CCD)  camera  technology. 
When  fully  calibrated,  sky  maps  of  structures  having 
enhanced  electron  density  in  the  inner  heliosphere  can  be 
produced.  Over  its  5-year  lifetime,  SMEI  has  observed  over 
300  ICMEs  traveling  through  the  inner  heliosphere.  Webb 
et  al.  [2006a]  summarized  the  SMEI  observations  of  ICMEs 
during  the  first  1.5  years  of  the  mission.  At  least  25  of  these 
transients,  typically  observed  as  frontside  halo  events  by  the 
SOHO  Large  Angle  Spectroscopic  Coronagraph  (LASCO), 
were  observed  by  SMEI  propagating  to  and  beyond  1  AU, 
and  were  associated  with  major  geomagnetic  storms  at  Earth. 


Copyright  2009  by  the  American  Geophysical  Union 
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Figure  1.  (top)  SMEI  all-sky  Aitoff  sky  map.  (bottom) 
A  single-orbit  difference  image  showing  both  a  ICME 
(shaded  arrows)  and  the  main  obscuration  effects 
viewed  by  SMEI  (open  arrows).  Lettered  arrows  are 
missing  data  resulting  from  A,  the  zone  of  exclusion 
near  Sun;  B,  three  camera  frames  that  are  shuttered 
because  of  bright  sunlight  (adjacent  white-black  areas 
have  sunlight-saturated  CCD  pixels  that  are  not  bright 
enough  to  close  the  shutter);  C,  particle  enhancements 
from  both  the  polar  zones  and  SAA;  and  D,  auroral 
light.  The  ICME  appears  as  arcs  seen  in  three 
unobscured  areas  of  the  inner  heliosphere.  This  was 
the  first  earthward-directed  halo  ICME  observed  by 
SMEI  on  29  May  2003  (see  Figure  2  and  3). 


[3]  The  successful  operation  of  SMEI  is  an  important  step 
toward  demonstrating  improved  space  weather  forecasts. 
Early  detection  of  potentially  damaging  Earth-directed 
ICMEs  can  help  protect  space  assets  and  maintain  stable 
communications.  A  geomagnetic  storm  initiates  a  wide 
variety  of  adverse  effects,  including  increases  in  trapped 
magnetospheric  particles,  degraded  satellite  communica¬ 
tion  and  surveillance  systems,  increased  drag  and  deteri¬ 
oration  of  satellite  altitude  control,  and  destructive  surges 
in  ground  power  grids.  For  example,  the  "Halloween" 
storms  in  October- November  2003  caused  everything 
from  power  grid  failures  to  satellite  shut  downs  (see  Webb 
and  Allen  [2004]  for  a  list  of  anomalies). 

[4]  SMEI's  mission  as  an  Air  Force  experiment  is  to  dem¬ 
onstrate  the  ability  to  detect  and  track  ICMEs  from  near  the 
Sun  to  Earth,  thus  providing  a  new  capability  for  forecast¬ 
ing  the  occurrence  of  storms.  SMEI  has  accomplished  this 
by  detecting  a  number  of  geoeffective  ICMEs  at  elongations 


of  20-30°  (from  Sun  center),  or  estimated  distances  as  far 
away  as  2/3  of  the  distance  between  the  Sun  and  Earth. 
Depending  on  the  speed  of  the  ICME  front,  these  distances 
would  correspond  to  lead  times  of  10  hours  to  1-2  days. 
SMEI  was  not  designed  to  be  an  operational  mission,  so 
most  SMEI  forecasts  were  retrospective,  with  a  couple  of 
exceptions  (examples  of  real-time  SMEI  forecasts  are  noted 
in  section  4.2.2).  But  these  results  do  prove  the  principle  that 
SMEI  could  detect  even  fast  Earthward  ICMEs  <1  day 
before  their  Earth  arrival  and,  thus,  that  a  similar  system, 
with  automated  event  detection  and/or  24/7  operator  mon- 


CDAW  STORM  NUMBER  63 


CDAW  STORM  NUMBER  63 


Figure  2.  The  29-30  May  2003  ICME  (Figure  1 
(bottom)).  CDAW  storm  63  (top)  elongation  versus 
time  and  (bottom)  distance  versus  time  plots.  See 
Figure  3.  The  symbols  are  as  follows:  squares,  leading 
edge  of  CMEs  in  LASCO  C2  and  C3;  triangles,  leading 
edge  of  ICMEs  in  SMEI;  asterisks,  shock  ATs  at  Ll  at 
ACE  or  Wind  (Shock  Date  and  Time  column  of  Table  2); 
crosses,  SMEI  aurora  onset  time  (Aurora  Date  column 
of  Table  2);  vertical  bars,  storm  onset  and  peak  times 
(Dst  Date  column  of  Table  2).  The  d-t  data  is  used  to 
calculate  the  arrival  times  of  the  ICME  at  1  AU  based 
on  linear  fits  to  the  three  separately  tracked  SMEI 
features;  the  two  later  ATs  are  given  in  Predicted  3-D 
column  of  Table  2. 
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ACHIEVE  L2  ACELEVEL2  Osi  KP  ICMES 


Figure  3.  ACE  solar  wind  and  geomagnetic  data  during  this  period,  from  29  to  31  May,  for  the  29- 
30  May  2003  ICME  (Figure  lb;  CDAW  storm  63).  The  plot  shows  (a)  the  magnetic  field  intensity, 
the  (b)  polar  and  (c)  azimuthal  IMF  components,  (d)  the  solar  wind  proton  temperature  (black)  and 
expected  proton  temperature  (red)  with  the  black  shading  indicating  where  the  proton 
temperature  falls  below  the  expected  temperature,  (e)  the  proton  density,  (f)  velocity,  (g)  helium/ 
proton  ratio,  (h)  the  observed  and  expected  07/06  ratio,  (i)  the  Bz  component  of  the  IMF  in  GSM 
coordinates,  the  (j)  Dst  and  (k)  Kp  geoindices,  and  (1)  a  0/1  value  indicating  whether  or  not  the 
signatures  suggest  an  ICME  flow  [e.g.,  Richardson  and  Cane ,  2004].  The  vertical  red  bar  indicates 
the  onset  time  of  bright  aurora  in  SMEI  (arrow)  when  the  ICME  front  was  at  an  elongation  of  80°. 
The  vertical  green  line  indicates  the  AT  of  the  first  shock  at  ACE.  ACE  plot  courtesy  of 
I.  Richardson. 


itoring,  and  adequate  data  latency  would  constitute  a 
reliable  tool  for  early  warning  of  storms. 

(5]  The  first  Earth-directed  ICME  detected  and  analyzed 
in  the  SMEI  data  occurred  in  late  May  2003  [Tap pin  et  ai, 
2004]  (Figures  1  -3).  This  fast  (projected  speed  of  1000  km/s) 
event  was  detected  at  ^30°  elongation,  about  15  hours 
before  it  passed  over  the  Earth  on  29-30  May  causing  a 


major  geomagnetic  storm.  SMEI  also  observed  ICMEs 
associated  with  a  series  of  large  solar  events  during  2  weeks 
in  late  October  and  early  November  2003.  These  led  to  high 
levels  of  energetic  particles  in  geospace  and  produced 
overlapping  large  geomagnetic  storms  on  29-30  October. 
The  28  and  29  October  flares  had  peak  X-ray  fluxes  among 
the  largest  ever  recorded,  shock  waves  among  the  fastest 
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ever  observed,  and  two  intense  storms  among  the  14 
strongest  storms  dating  back  to  1936.  SMEI  observed  the 
first  of  these  halo  ICME  events,  again  starting  at  distances 
about  1/3  of  the  way  from  the  Sun  to  Earth,  10  hours  before 
Earth  arrival. 

[6]  The  focus  of  this  paper  is  on  comparison  of  the  SMEI 
ICME  observations  with  the  intense  geomagnetic  storms 
that  were  the  focus  of  the  NASA  Living  With  a  Star  (LWS) 
Geostorm  Coordinated  Data  Analysis  Workshop  (CDAW). 
This  CDAW  brought  together  scientists  from  many  LWS 
disciplines  to  interact  and  collaborate  on  questions  related 
to  CMEs,  ICMEs,  and  their  connections  to  geomagnetic 
and  ionospheric  storms.  The  focal  point  of  the  CDAW  was 
the  set  of  all  major  geomagnetic  storms  (defined  as  having 
a  peak  Dst  <  -100  nT)  of  cycle  23,  during  the  SOHO  era 
from  1996  to  2004.  Two  Geostorm  CDAWs  were  held,  in 
2005  and  2007;  the  purpose  of  the  second  was  to  extend  the 
list  of  storms  through  2005  and  to  have  a  scientific  meeting 
to  present  first  results.  The  data  and  some  preliminary 
results  from  the  CDAWs  are  available  online  at  http:// 
cdaw.gsfc.nasa.gov/geomag  cdaw /,  especially  under 
"Data"  and  "Resources." 

(7]  There  were  four  CDAW  working  groups;  working 
group  1  (WG1)  on  Drivers  of  Geomagnetic  Storms,  WG2 
on  Geomagnetic  Storm  Mechanisms,  WG3  on  Ionospheric 
Storms,  and  WG4  on  Prediction  of  Geomagnetic  Storms. 
The  main  goal  of  WG1  was  to  try  to  identify  and  categorize 
the  solar  and  interplanetary  sources  of  each  of  the  intense 
storms.  The  major  results  of  this  process  were  presented 
by  Zhang  et  al.  [2007a,  2007b].  A  second  paper  by  Richardson 
etal.  [2006,  2007]  summarized  results  for  the  ten  storms  that 
were  associated  with  corotating  interaction  regions  (CIRs). 
Since  initiating  operations  in  February  2003,  SMEI  was  able 
to  observe  the  associated  ICMEs  for  89%  of  the  intense 
storms.  Here  we  describe  the  SMEI  observations  and  anal¬ 
yses  for  the  26  CDAW  storms  that  occurred  during  the 
SMEI  observations  from  May  2003  to  September  2005,  as 
well  as  the  two  other  intense  storms  that  occurred  after  that 
period  through  the  end  of  2007. 

[fi]  In  section  2  we  summarize  the  SMEI  mission,  empha¬ 
sizing  its  observations  of  ICMEs  that  were  likely  directed 
toward  the  vicinity  of  the  Earth.  In  section  3  we  describe  the 
main  results  comparing  the  SMEI  ICME  observations  with 
the  intense  storms,  and  in  section  4  the  implications  of  the 
SMEI  results  for  the  forecasting  of  such  storms  at  Earth.  The 
results  are  summarized  and  discussed  in  section  5. 

2.  SMEI  Mission  and  How  We  Identify 

Earthward  ICMEs 

2.1.  SMEI  Background 

[9]  SMEI  was  launched  in  January  2003  on  the  Air  Force 
Space  Test  Program's  Coriolis  satellite.  Coriolis  is  in  a 
dawn-dusk.  Sun-synchronous,  circular  polar  orbit  along 
the  Earth's  terminator  at  an  altitude  of  840  km  and  an 
inclination  of  98°  relative  to  the  equatorial  plane,  and 
maintains  a  zenith-nadir  orientation.  SMEI  has  three  baf¬ 


fled  cameras,  each  covering  a  narrow  3°  x  60°  strip  of  the 
sky.  Camera  1  points  to  the  night  sky.  Camera  2  is  centered 
on  the  terminator  and  Camera  3  points  near  the  Sun.  The 
cameras  are  mounted  on  the  satellite  with  their  fields  of 
view  (FOVr)  aligned  end  to  end  and  slightly  overlapping,  so 
that  the  FOV  swath  is  a  3°  wide  strip  extending  160°  along 
an  approximate  great  circle  with  the  ends  near  the  orbit 
axis.  Each  camera  takes  continuous  4-s  exposures  in  broad¬ 
band  white  light  during  each  orbit  with  14  orbits  per  day. 
Since  the  cameras  are  fixed  to  the  spacecraft,  they  image 
nearly  90%  of  the  entire  sky  during  each  orbit.  Gaps  in 
coverage  include  a  zone  of  exclusion  of  ^20°  radius  about 
the  sunward  orbital  pole,  a  smaller  circle  in  the  opposite 
direction,  and  occasional  areas  shuttered  because  of  sun¬ 
light  in  the  sunward  camera.  The  instrument  operates 
continuously,  so  the  primary  data  product  comprises  a 
sequence  of  14  all-sky  images  per  day.  It  has  maintained  a 
duty  cycle  of  85%,  interrupted  only  for  periodic  calibration 
and  diagnostic  purposes  and  occasional  software  anomalies 
and  telemetry  problems.  SMEI  views  the  outward  flow  of 
density  structures  emanating  from  the  Sun  by  observing 
Thomson-scattered  sunlight  from  heliospheric  plasma.  The 
primary  objective  of  SMEI  is  to  demonstrate  the  feasibility 
of  using  such  instrumentation  to  forecast  the  arrival  of  the 
ICMEs  at  Earth.  SMEI  has  achieved  this  goal  by  imaging 
nearly  the  entire  sky  every  102  min,  the  spacecraft's  orbital 
period,  with  sufficient  sensitivity  and  photometric  stability 
to  be  able  to  detect  faint  transient  disturbances  against  the 
much  brighter,  but  relatively  unchanging  stellar  and  zodi¬ 
acal  backgrounds. 

[10]  The  individual  3°  x  60°  data  frames  are  combined 
on  the  ground  to  form  composite  heliospheric  all-sky  maps. 
We  typically  use  two  types  of  Sun-centered  projection:  an 
equal-area  Hammer-Aitoff  projection  and  an  angular 
"fisheye"  projection  wherein  the  distance  from  the  center 
of  the  image  is  proportional  to  the  angle  from  the  camera 
view  direction.  An  Aitoff  sky  map  provides  a  reasonably 
undistorted  representation  of  the  sunward  hemisphere, 
where  most  heliospheric  disturbances  are  best  observed, 
but  at  the  cost  of  more  serious  distortion  around  the 
perimeter  of  the  image  region.  Figure  1  (top)  shows  an 
Aitoff  sky  map,  with  ecliptic  north  at  the  top  and  east  to 
the  left.  Fisheye  projections  (Figure  5  shows  an  example) 
have  relatively  little  distortion  out  to  ~100°  elongation  from 
the  Sun  and  are  particularly  effective  in  following  Earth- 
directed  events. 

[11]  Details  of  the  SMEI  data  handling,  processing, 
background  light  subtraction  and  sky  map  production, 
particularly  their  pertinence  for  observing  ICMEs,  are 
discussed  by  Jackson  et  al.  [2004]  and  Webb  et  al.  [2006a]. 
Basically,  SMEI  is  designed  to  measure  heliospheric 
Thomson-scattered  light  integrated  along  each  line  of 
sight,  like  any  white  light  coronagraph.  Since  the  bright¬ 
ness  of  this  emission  falls  off  rapidly  with  distance  from  the 
Sun,  care  is  needed  in  generating  the  full-sky  SMEI  images 
to  insure  that  the  variations  in  brightness  from  these  very 
faint  signals  are  detectable  over  the  entire  sky. 
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[12]  There  are  a  number  of  artifacts  or  backgrounds  that 
can  affect  the  SMEI  sky  maps.  When  a  SMEI  camera  points 
close  to  the  Sun,  sunlight  illuminates  inner  portions  of  the 
baffle  and  can  saturate  the  CCD  exposure.  The  sunward 
edge  of  the  Camera  3  FOV  is  offset  20°  from  the  sunward 
pole  of  the  orbital  plane  creating  a  circular  zone  of  avoid¬ 
ance  ^20°  in  radius  on  the  orbit  maps.  Because  Coriolis  is  in 
a  Sun-synchronous  orbit  and  not  Sun  pointed,  this  exclu¬ 
sion  zone  moves  north  and  south  of  Sun  center  during  the 
year.  To  protect  the  CCD,  a  shutter  in  the  light  path  closes 
when  the  solar  illumination  is  too  high.  The  sunward 
Camera  3  operates  at  much  higher  temperatures  than 
intended,  lowering  its  overall  sensitivity  primarily  because 
of  particle-induced  radiation  damage  on  the  CCD.  The 
background  from  individual  stars,  the  extended  emission 
from  the  Galactic  plane,  and  the  zodiacal  light  are  each  far 
brighter  than  the  faint  heliospheric  structures  (ICMEs)  that 
SMEI  is  designed  to  detect.  Other  artifacts  include  the 
residuals  from  bright  stars,  data  dropouts,  and  saturated 
regions  due  to  the  Moon  and  bright  planets.  Another 
artifact  that  is  caused  by  the  interactions  of  trapped  radia¬ 
tion  belt  particles  with  the  CCD  and  its  processing  elec¬ 
tronics  occurs  on  nearly  every  orbit  and  can  affect  large 
numbers  of  CCD  pixels.  These  particles  typically  appear  as 
the  satellite  crosses  the  polar  auroral  zones  and  during 
passages  through  the  South  Atlantic  Anomaly  (SAA).  On 
the  maps  these  regions  appear  as  large  area  bands  that  have 
a  granular  appearance. 

[13]  Another  type  of  "background"  is  an  occasional 
visible  light  phenomenon  associated  with  the  geoaurora 
[Mizuno  et  al,  2005].  This  phenomenon  occurs  in  the  auroral 
zones  and  polar  regions  above  Coriolis's  orbit  (>840  km); 
SMEI's  observations  of  this  high-altitude  aurora  were  un¬ 
anticipated  and  constitute  a  major  discovery.  This  auroral 
light  viewed  by  SMEI  can  be  very  bright,  and  it  is  well 
correlated  in  both  brightness  and  frequency  with  periods  of 
enhanced  geomagnetic  activity  when  the  indices  Kp  and  Dst 
have  peak  values  >^6  and  <  -60  nT,  respectively.  As  such  it 
has  proven  to  be  of  great  value  from  a  space  weather 
perspective  in  identifying  the  time  of  arrival  of  geoeffective 
shocks  and  ejecta  and  the  onset  and  duration  of  the  asso¬ 
ciated  storms. 

[14]  Figure  1  (bottom)  illustrates  these  artifacts  (except 
for  the  moon)  along  with  a  true  ICME  on  an  orbital  differ¬ 
ence  image.  The  open  lettered  arrows  point  to  the  artifacts. 
The  ICME  (shaded  arrows)  appears  as  arcs  in  three  unob¬ 
scured  areas  of  the  inner  heliosphere.  This  was  the  first 
Earthward-directed  halo  ICME  observed  by  SMEI,  on 
29  May  2003. 

2.2.  Identifying  Earthward  ICMEs  in  SMEI 

[is]  The  primary  backgrounds,  which  are  from  individ¬ 
ual  stars,  the  Galactic  plane  emission  and  the  zodiacal  light, 
can  be  removed  to  first  order  by  producing  orbital  differ¬ 
ence  images,  so-called  "running  differences."  In  effect,  a 
constant  heliospheric  signal  cancels  out  in  the  difference 
map,  leaving  only  the  changes  in  heliospheric  structures 


from  one  orbit  to  the  next.  We  found  that  ICMEs  are  usually 
too  faint  to  be  detectable  on  the  direct  maps  or  movies  made 
from  them,  but  become  readily  observable  as  moving 
structures  on  the  difference  movies.  More  advanced  tech¬ 
niques  are  being  developed  at  UCSD  and  AFRL  to  remove 
these  background  signals  and  better  calibrate  the  SMEI 
data  [e.g.,  Jackson  et  al.,  2006;  Hick  et  al,  2005],  Such  pro¬ 
cessed  sky  maps  have  now  been  produced  for  most  of  the 
SMEI  mission  and  are  being  used  in  the  data  analysis  on  a 
case-by-case  basis.  We  have  used  these  images  and  movies 
to  analyze  the  storm  events  herein. 

[ie]  Webb  et  al  [2006a]  analyzed  the  SMEI  observations 
of  ICMEs  for  the  first  1.5  years  of  the  SMEI  mission.  The 
SMEI  ICME  occurrence  rate  was  about  1/3  CMEs/day, 
about  an  order  of  magnitude  less  than  the  LASCO  CME 
rate  for  the  same  time  period.  Thus,  SMEI  sees  only  a 
fraction  of  all  eruptions  identified  as  CMEs  in  the  LASCO 
data,  and  we  are  currently  comparing  CMEs/ICMEs  ob¬ 
served  with  both  instruments  to  characterize  the  differ¬ 
ences.  Some  reasons  for  the  lower  SMEI  rate  are  that  SMEI 
is  less  sensitive  to  CMEs  erupting  off  the  backside  of  the 
Sun  because  of  Thomson  scattering  at  these  large  elonga¬ 
tions,  and  frequent,  multiple  CMEs  can  merge  at  large 
distances  from  the  Sun.  Recently,  Howard  and  Simnett  [2008] 
examined  appropriate  LASCO  data  for  evidence  of  CMEs 
corresponding  to  189  SMEI  events  observed  from  February 
2003  to  September  2005.  Of  those  SMEI  ICMEs,  143  (76%) 
were  associated  with  one  or  more  LASCO  CMEs,  31  (16%) 
were  only  weakly  associated  with  LASCO  activity,  and  15 
(8%)  had  no  apparent  LASCO  activity.  That  SMEI  observes 
ICME  material  in  the  inner  heliosphere  that  is  not  visible 
near  the  Sun  is  an  important  finding  that  is  currently  under 
study. 

[17]  Webb  et  al  [2006a]  reported  that  the  average  dura¬ 
tion  of  SMEI  ICMEs  was  ~16  hours  with  a  maximum  of 
about  three  days.  The  SMEI  ICME  angular  spans  or  widths 
averaged  42°,  ranging  up  to  107°,  but  this  is  a  lower  limit 
because  of  obscuration  of  the  sides  of  many  ICMEs.  The 
average  peak  brightness  of  a  SMEI  ICME  was  2.3  S10  units 
(an  S10  unit  is  the  equivalent  flux  of  10th  magnitude, 
solar-type  stars  per  square  degree).  The  maximum  ob¬ 
served  elongation  to  which  SMEI  could  track  ICMEs  was 
^145°.  The  SMEI  ICMEs  could  be  grouped  into  seven 
interpretative  categories,  of  which  the  most  common  were 
limb  ICMEs,  distant  wide  arcs  and  Earthward-directed 
ICMEs.  The  latter  are  of  the  most  interest  for  this  paper. 

[is]  Figures  4a -4d  are  an  example  of  the  limb  ICME 
class,  one  of  two  similar,  fast,  west  limb  loop  events  likely 
associated  with  NOAA  active  region  0365,  the  region 
associated  with  the  earlier  28-29  May  2003  halo  ICME 
(Figures  1-3)  when  the  active  region  was  near  Sun  center. 
The  limb  event  on  31  May  led  to  a  circular  loop  CME  seen 
in  LASCO  C2  and  C3  images  (inset)  and  the  SMEI  ICME. 
The  loop's  location  and  width  (^60°)  were  similar  in  both 
the  LASCO  and  SMEI  views.  Figure  4e  is  the  elongation 
versus  time  plot  of  the  ICME  front.  Assuming  projection 
in  the  skyplane,  the  speeds  of  the  front  calculated  from  the 
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Figure  4.  (a-d)  Example  of  a  loop  CME  viewed  off  the  solar  limb  by  LASCO  (blue)  and  SMEI 
Camera  3  (red).  The  event  on  31  May  2003  was  associated  with  a  2B,  M9.3  flare  with  onset  at 
0215  UT  and  disk  location  S07°W65°.  This  led  to  a  circular  loop  CME  seen  in  LASCO  (C3  image  in 
inset),  (e)  Elongation  versus  time  plot  of  the  ICME  front  seen  in  SMEI  with  27°  <  5  <  44°.  The  linear 
fit  extrapolation  earlier  in  time  matches  well  with  the  flare  (F)  onset  and  LASCO  C2,  C3  height 
times.  From  Webb  et  al.  [2006a]. 


LASCO  (1830  km/s)  and  SMEI  (^1450  km/s)  data  agree 
reasonably  well.  These  events  illustrate  that  at  least  some 
ICMEs  appear  to  undergo  little  change  as  they  propagate 
outward  from  their  low  coronal  origins,  in  this  case  out  to 
45°  elongation.  Such  bright,  relatively  fast  ICMEs  occur¬ 
ring  near  the  solar  limb  are  fairly  easy  to  detect  in  SMEI's 
sunward  camera  and  to  correlate  with  solar  surface  events 
and  LASCO  CMEs.  This  is  because  the  intensity  of  light 
emitted  via  the  Thomson  scattering  process  is  maximized 
for  plasma  that  is  closest  to  the  Sun  along  the  line  of  sight. 
However,  in  the  SMEI  images  such  limb  CMEs  fade 
rapidly  as  they  move  outward,  as  expected  for  Thomson 
scattering  |see,  e.g.,  Vourlidas  and  Howard ,  2006).  Likewise, 
this  effect  implies  that  SMEI  will  not  detect  much,  if  any 
emission  from  events  ejected  from  the  backside  of  the  Sun. 

[is|  The  distant  arc  events  are  fairly  common,  though 
often  faint  and  subtle  in  individual  frames,  and  comprise 
about  1/3  of  all  the  SMEI  events.  Because  of  the  Thomson- 


scattering  geometry,  these  distant  arcs  could  either  be 
directed  toward  the  Earthward  hemisphere  of  the  sky  or 
in  some  cases  could  be  near  Earth.  Typical  examples  are 
the  events  in  Figures  5  and  6.  The  arc  shown  in  Figure  5a 
swept  beyond  90°  elongation  through  Camera  2,  which  is 
centered  on  the  terminator,  into  Camera  1,  thus  from  the 
dayside  into  the  night  sky,  at  about  the  time  a  geomag¬ 
netic  storm  began.  Such  terminator-crossing  events  are 
best  observed  in  the  southeast  quadrant  of  the  SMEI  FOV, 
where  the  obscuration  from  particles  in  the  auroral  zones 
and  the  SAA  are  minimal. 

[ 2o[  Since  ICMEs  are  the  primary  cause  of  major  space 
weather  events,  such  as  storms  and  particle  enhancements, 
it  is  important  to  be  able  to  detect  and  track  their  fronts  as 
they  propagate  through  the  inner  heliosphere  toward 
Earth.  Until  the  launch  of  SMEI,  forecasters  had  almost 
no  information  on  the  movement  of  a  geoeffective  ICME 
once  it  left  the  FOV  of  a  coronagraph,  such  as  the  LASCO 
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Figure  6.  The  3-10  November  2004  LASCO  CMEs, 
SMEI  ICMEs  (Figure  5),  and  time  markers  at  1  AU  for 
CDAW  storms  77  and  78.  (top)  Elongation-time  plot 
and  (bottom)  distance-time  plot  Symbols  are  the  same 
as  in  Figure  2. 

C3  which  views  to  a  radius  of  7.5°,  or  30  R$,  from  the  Sun. 
Forecasters,  however,  have  been  aided  by  LASCO's  ability 
to  frequently  observe  halo  CMEs.  Such  CMEs  appear  as  a 
bright  ring  surrounding  the  occulter  indicating  their  direc¬ 


tion  along  the  Sun-Earth  line;  additional  observations  of 
frontside  surface  activity  are  used  to  indicate  whether  it  is 
Earth  directed.  Several  radio  techniques  have  been  used  to 
remotely  detect  and  track  disturbances  in  the  interplane¬ 
tary  medium.  These  are  kilometric  radio  observations  from 
space  and  interplanetary  scintillation  (IPS)  observations 
from  the  ground.  The  kilometric  observations  can  track 
the  emission  typically  only  from  strong  shocks  traveling 
ahead  of  fast  ICMEs.  The  IPS  technique  suffers  from 
relatively  poor  temporal  (^24  hours)  and  spatial  resolution, 
ionospheric  noise,  and  a  model  dependence  for  interpret¬ 
ing  the  signal  as  density  or  mass. 

(21]  SMEI  has  now  observed  at  least  30  ICMEs  that  were 
likely  directed  toward  the  general  vicinity  of  Earth.  These 
are  typically  a  subset  of  the  arc-like  category  of  events. 
Those  most  likely  Earth  directed  have  these  characteristics: 
(1)  They  are  associated  with  solar  activity  near  Sun  center 
and  with  LASCO  halo  CMEs;  (2)  they  tend  to  be  arc-like  and 
wide,  with  spans  on  the  order  of  60°;  (3)  the  arcs  appear  to 
sweep  by  Earth's  terminator  from  the  dayside  into  the  night 
side,  i.e.,  they  move  across  the  Camera  2  FOV;  and  (4)  the 
timing  of  this  terminator  passage  tends  to  correspond  to  the 
onset  of  auroral  brightening  in  SMEI.  Examples  of  such 
events  are  shown  in  Figures  1,  5,  and  11. 

3.  Review  of  Prior  Studies  Comparing  SMEI 
ICMEs  With  Geomagnetic  Storms 

[22J  There  have  been  several  previous  studies  compar¬ 
ing  SMEI  ICMEs  with  geomagnetic  storms  that  are  perti¬ 
nent  to  this  paper.  Since  not  all  of  these  studies  have  been 
published,  we  summarize  their  results  in  this  section.  The 
studies  are  grouped  into  three  topical  areas.  The  common 
study  among  all  three  is  Howard  et  al.'s  [2006]  study  of 
20  SMEI  ICMEs  compared  with  shocks  at  ACE  at  LI.  That 
study  was  basically  the  precursor  to  this  paper. 

3.1.  Association  of  SMEI  ICMEs  With  Moderate 
to  Major  Storms 

[23]  Two  preliminary  and  independent  studies  were 
made  of  SMEI's  capability  to  detect  and  track  ICMEs  that 
subsequently  caused  major  storms  at  Earth.  These  results 


Figure  5.  SMEI  ICMEs  (CDAW  storms  77  and  78).  Four  views  of  SMEI  ICMEs  in  orbit  difference  images  ((a)  5 
November  2004  at  1049  UT,  first  loop  appears  to  east  (arrow)  in  Cameras  3  and  2;  (b)  7  November  2004  at  0832  UT, 
two  arcs  near  Earth  (arrows);  (c)  7  November  2004  at  1337  UT,  onset  of  aurora  (large  black  and  white  areas  in  upper 
left  and  lower  right);  and  (d)  9  November  2004  at  0110  UT,  loop  to  northwest  (arrow)  in  Camera  2).  The  plus  sign 
marks  the  Sun's  location.  Note  that  the  SMEI  FOV  is  shifted  northward  during  this  season.  ACE  solar  wind  data 
plot  shows  (e)  the  magnetic  field  intensity,  the  (0  polar  and  (g)  azimuthal  IMF  components,  (h)  the  solar  wind 
proton  temperature  (black)  and  expected  proton  temperature  (red)  with  the  black  shading  indicating  where  the 
proton  temperature  falls  below  the  expected  temperature,  (i)  the  proton  density,  (j)  velocity,  (k)  helium/proton 
ratio,  (1)  the  observed  and  expected  07/06  ratio,  (m)  the  Bz  component  of  the  IMF  in  GSM  coordinates,  the  (n)  Dst 
and  (o)  Kp  geoindices,  and  (p)  a  0/1  value  indicating  whether  or  not  the  signatures  suggest  an  ICME  flow  [e.g., 
Richardson  and  Cane,  2004J.  The  vertical  red  bar  indicates  the  onset  time  of  bright  aurora  in  SMEI  (arrow)  on  7  and 
9  November  2004  when  the  ICME  fronts  were  at  elongations  of  105  and  85°,  respectively.  The  second  ICME  may  be 
more  closely  associated  with  the  first  Dst  dip  (D)  before  the  storm  78  peak.  The  vertical  green  lines  indicate  the  ATs 
of  shocks  at  ACE  (see  text).  ACE  plot  courtesy  of  1.  Richardson. 
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Figure  7.  Example  of  the  comparison  of  ICMEs  and  ACE  shocks  from  Howard  et  al.'s  [2006]  study. 
Event  of  21-23  January  2004  (CDAVV  storm  73).  (a)  LASCO  C3  halo  CME,  (b)  four  views  of  the 
SMEI  transient,  and  (c)  SMEI  elongation  versus  time  plot  with  examples  of  various  fits.  These 
allowed  Howard  et  al.  [2006]  to  predict  the  ICME  arrival  time  and  compare  it  with  the  actual  (shock) 
arrival  time.  The  arrow  indicates  the  onset  time  of  SMEI  aurora. 


were  first  reported  by  Johnston  et  al.  [2004]  and  discussed 
by  Howard  et  al.  [2006]  and  Webb  et  al.  [2006a].  In  one  study 
we  included  all  moderate  level  or  greater  storms,  defined 
as  having  a  peak  Dst  <  -60  nT,  over  an  approximately  2-year 
time  period.  For  85%  (39  of  46)  of  these  storms,  SMEI 
detected  a  ICME  within  2  days  prior  to  the  storm  onset. 
Whether  all  the  storms  were  ICME  driven  or  whether  the 
observed  ICMEs  could  be  clearly  associated  with  the  sub¬ 
sequent  storms  was  not  established,  nor  were  the  charac¬ 
teristics  of  those  ICMEs  examined. 

[24]  In  another  study  we  examined  the  sources  of  only 
the  most  intense  (peak  Dst  <  -100  nT)  geomagnetic 
storms  during  the  first  2  years  of  SMEI  observations,  from 
February  2003  to  2005.  There  were  14  such  storms  and 
SMEI  had  suitable  data  during  12  of  them.  For  10  of  those 
12  storms  (83%),  SMEI  observed  associated  Earthward- 
directed  ICMEs.  During  all  12  storms  SMEI  also  observed 
the  bright  auroral  light  associated  with  the  storm.  The 
mean  time  difference  between  the  first  SMEI  observation 
of  the  ICxME  and  the  arrival  of  the  associated  shock  at  Earth 
was  18.6  hours.  The  mean  difference  between  the  first 
ICME  observation  and  the  storm  onset  was  29.25  hours. 


[25]  These  two  studies  were  the  first  to  show  that  SMEI 
could  detect  fairly  fast  Earthward  ICMEs  from  ^10  hours 
to  ~2  days  before  their  Earth  arrival  and  that,  given  much 
better  data  latency,  a  future  SMEI-type  heliospheric  im¬ 
ager  could  be  used  to  forecast  the  onset  and  maybe  even  the 
intensity  of  geomagnetic  storms.  These  preliminary  studies 
have  now  been  expanded  and  combined  with  the  identifi¬ 
cation  results  of  the  Geostorm  CDAW  to  cover  all  intense 
storms  during  the  SMEI  observations  from  February'  2003  to 
2007.  Those  results  are  the  focus  of  this  paper. 

3.2.  Comparing  LASCO  CMEs  and  SMEI  ICMEs 
With  ACE  Shocks 

[26]  Howard  et  al.  [2006]  studied  20  ICMEs  observed  by 
SMEI  over  its  first  19  months  that  were  associated  with 
LASCO  halo  CMEs  and  shocks  observed  by  ACE  at  the 
LI  point.  Their  study  focused  on  computing  distance-time 
profiles  using  both  the  LASCO  CME  and  SMEI  ICME  data 
to  determine  ICME  speed,  and  then  to  predict  its  time  of 
arrival  at  Earth.  They  used  the  shock  arrival  time  (AT)  at  LI 
as  the  indicator  for  ICME  arrival,  thus  ignoring  the  shock 
standoff  distance  (time)  from  the  ejecta.  They  then  com- 
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Figure  8.  Time  difference  histogram  comparing  the  predicted  (LASCO  and  SMEI  Ejecta)  and 
actual  (ACE  shock)  arrival  times  (15  events  with  storms).  The  colors  denote  different  levels  of 
storms,  i.e.,  green  are  small,  yellow  are  medium,  and  red  are  large  storms,  respectively.  The  cross- 
hatched  bars  are  for  storms  with  no  sudden  storm  commencements  (SSCs);  see  Howard  et  al  [2006] 
for  details.  Algebraic  mean  is  +4.0  hours.  RMS  is  11.4  hours. 


pared  the  predicted  ICME  AT  with  the  actual  shock  ATs  at 
ACE  at  LI. 

[27]  The  angular  distance  from  the  Sun  traversed  by  an 
ICME  is  typically  characterized  by  its  solar  elongation 
angle,  t.  In  ecliptic  coordinates,  the  center  of  the  Sun  is 
at  s  =  0°,  e  is  the  polar  distance  and  position  angle  (PA)  is 
the  associated  polar  coordinate.  For  SMEI,  the  elonga¬ 
tion  5  of  the  leading  edge  of  an  ICME  versus  time  plot 
provides  the  primary  information  on  the  outward 
motion  of  the  ICME.  Only  angular  speeds  can  be  de¬ 
rived  from  fits  to  the  e  versus  time  plot  for  a  given  ICME. 
To  calculate  true  speeds,  some  assumptions  or  a  partic¬ 
ular  model  must  be  used.  Howard  et  al.  [2006]  computed 
speeds  on  the  basis  of  the  so-called  "point  P  approxi¬ 
mation."  This  depends  on  the  assumption  that  the  ICME 
front  can  be  approximated  as  a  Sun-centered  circular  arc 
and  that  the  maximum  density  along  the  line  of  sight 
occurs  at  the  closest  approach  point  of  that  line  of  sight 
to  the  Sun,  i.e.,  at  the  tangent  point,  P,  of  the  arc.  The 
distance  from  the  Sun  to  P  then  is  equal  to  the  sine  of  the 
observed  5,  since  SMEI  is  at  1  AU  from  the  Sun,  and 
the  speed  is  the  change  of  this  distance  with  time, 
converted  to  appropriate  units,  such  as  km /sec.  The 
point  P  approximation  provides  only  lower  limits  to 
the  true  distances  and  speeds  of  ICMEs  as  viewed  by 
SMEI  and  is  less  accurate  for  ICMEs  approaching  Earth. 

[28]  Figure  7  shows  an  example  of  one  of  the  20  Howard  et 
al.  [2006]  events,  their  event  number  15  on  21  -22  January 
2004.  YVe  feature  it  here  because  it  is  also  one  of  the  CDAYV 
storm  events,  number  73,  and  illustrates  the  method  they 
used.  Figure  7a  shows  the  LASCO  halo  CME  on  21  January, 
and  Figure  7b  shows  the  ICME  front  traveling  through  the 
SMEI  field  on  four  sky  maps  on  the  next  day.  Figure  7c 
shows  the  elongation-time  plot  with  the  LASCO  CME  and 
SMEI  ICME  points  and  the  ACE  shock  arrival  time  at  1  AU. 
A  linear  fit  through  the  SMEI  points  on  the  (point  P) 


distance-time  plot  (not  shown)  projects  to  2106  UT  on 
22  January,  whereas  the  actual  shock  arrival  time  was  at 
0105  UT  on  22  January.  So,  in  this  case,  SMEI  would  have 
predicted  a  storm  (shock)  time  ~20  hours  too  late. 

[2s]  During  2003-2004,  20  shocks  were  associated  with 
SMEI  ICMEs  with  15  of  these  also  associated  with  moder¬ 
ate  to  strong  storms.  Using  the  above  technique,  Howard  et 
al.  [2006,  Table  3]  determined  the  time  differences  between 
the  predicted  (LASCO  CMEs  and  SMEI  ICMEs)  and  actual 
(ACE  shocks)  ATs  for  these  15  events.  Figure  8  is  a  histo¬ 
gram  of  these  time  differences  between  the  predicted  and 
the  actual  ATs.  The  algebraic  mean  for  the  ICME  AT-shock 
time  differences  that  we  calculated  from  these  results  was 
+4.0  hours.  In  section  4.2.1  we  will  compare  these  results 
with  those  from  our  analysis  of  the  CDAYV  events. 

3.3.  Space  Weather  Forecast  Model  Corrections 
Using  SMEI 

[vo]  SMEI  provides  important  observations  of  the  prop¬ 
agation  and  evolution  of  ICMEs  which  contribute  to  under¬ 
standing  ICMEs  and  their  interaction  with  the  background 
solar  wind.  Timely  SMEI  observations  of  ICMEs  can  im¬ 
prove  space  situational  awareness  and  solar  wind  forecast 
skill.  At  a  fundamental  level,  SMEI  images  can  confirm  or 
refute  whether  an  ICME  that  is  predicted  by  a  physics- 
based  solar  wind  model  is  actually  en  route  to  Earth.  SMEI 
data  can  also  enable  midcourse  updates  or  corrections  to 
operational  solar  wind  forecasts  more  than  a  day  or  two 
before  the  disturbance  hits  Earth. 

[n]  Currently,  we  are  investigating  the  use  of  SMEI 
observations  to  compare  with  the  Hakamada-Akasofu- 
Fry  Version  2  (HAFv2)  3-D  solar  wind  model  [Hakamada 
and  Akasofu,  1982;  Fry  et  al,  2001,  2003].  The  HAFv2  model 
is  a  kinematic  model  that  predicts  solar  wind  conditions 
(speed,  density,  dynamic  pressure,  and  interplanetary 
magnetic  field)  at  Earth  and  elsewhere  in  the  solar  system 
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Figure  9.  (a  and  b)  Comparison  of  shock/compressed  IP  density  of  HAFv2  model  and  SMEI 
observations  for  the  29  May  2003  halo  ICME  (Figure  1),  associated  with  CDAW  storm  63;  (c)  HAFv2 
forecast  of  ecliptic  plane  IMF  showing  the  shock  propagating  through  a  preexisting  corotating 
interaction  region  (CIR).  The  reconstructed  HAF  density  maps  suggest  that  the  transient  viewed  by 
SMEI  is  brighter  to  the  east  because  of  the  added  contribution  from  the  CIR.  After  Sun  et  al.  [2008j. 


based  on  observations  at  the  Sun.  These  solar  wind 
parameters  are  key  inputs  to  numerical  prediction  models 
for  forecasting  space  weather  disturbances  that  impact  space 
operations.  The  IiAFv2  model  is  a  single-fluid,  Parker-like, 
modified  kinematic  model  that  projects  fluid  parcels  out¬ 
ward  from  the  rotating  Sun  along  fixed  radials  at  successive 
time  steps,  in  an  inertial  frame.  The  parcels  move  outward 
with  different  speeds  along  fixed  radials  as  the  Sun  rotates 
beneath  the  radial  grid.  The  frozen  field  condition  man¬ 
dates  that,  along  each  radial,  the  faster  parcels  do  not  move 
through  or  pass  slower  parcels.  Therefore,  the  kinematic 
flow  is  modified,  in  that  the  fluid  parcel  positions  are 
adjusted  to  account  for  the  stream-stream  interaction  as 
fast  parcels  (solar  wind  streams)  overtake  slower  ones.  The 
model  has  been  calibrated  so  that  HAFv2  output  matches  a 
1-D  MHD  solution.  If  the  speed  gradient  along  a  radial  is 
steep  enough,  corotating  interactions  and  interplanetary 
shocks  are  formed. 

[32]  The  HAFv2  model  uses  two  sets  of  inputs.  The  first 
set  establishes  the  steady  state  boundary  conditions  for 
the  “quiet,"  or  background  solar  wind,  and  the  second  set 
determines  the  time-dependent  boundary  conditions  for 
the  event-driven  solar  wind.  The  HAFv2  model  uses 
source  surface  maps  derived  from  synoptic  solar  observa¬ 


tions  [Arge  and  Pizzo,  2000;  Arge  et  al. ,  2004]  to  establish  the 
background  solar  wind  conditions.  It  is  also  driven  by 
proxies  for  energy  released  during  solar  events  to  model 
CMEs  and  interplanetary  shock  propagation.  The  HAFv2 
model  maps  the  disturbed  and  the  undisturbed  solar 
wind,  so  it  is  applicable  to  all  phases  of  the  solar  cycle. 
In  addition,  HAFv2  produces  chronological  sequences  of 
the  ecliptic  plane  interplanetary  magnetic  field  (IMF)  and 
other  solar  wind  parameters. 

(33)  We  are  using  HAFv2  to  compute  time-dependent 
solar  wind  density  maps  at  times  corresponding  to  SMEI 
image  frames.  We  compute  synthetic  skymaps  of  Thomson- 
scattering  brightness  by  integrating  heliospheric  density 
along  lines  of  sight  centered  at  Earth,  with  the  Sun  at  the 
center.  Total  brightness  was  computed  using  algorithms 
provided  by  one  of  us  (BVJ  [see  also  Jackson  et  al.,  2006; 
Billings,  1966]),  using  a  technique  developed  by  Sun  et  al. 
[2008J.  The  result  is  a  series  of  sky  plane,  Aitoff-projected 
images  of  the  modeled  ICME  taken  at  the  appropriate 
time.  These  were  converted  into  files  in  the  same  FITS 
format  as  those  from  SMEI,  so  that  measurements  of 
predicted  ICME  morphology  and  propagation  can  be 
obtained  in  the  same  manner  using  the  SMEI  analysis 
software. 
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[34]  Specific  studies  have  been  done  comparing  SMEI 
and  HAFv2  for  events  in  February  2003  and  February  2004 
l Howard  et  al ,  2007],  and  for  the  28-29  May  2003  event 
[Sun  et  al.,  2008].  Figure  9  illustrates  this  comparison  for 
the  latter  event.  Jackson  et  al  [2008]  describe  a  3-D  reconstruc¬ 
tion  of  this  event  using  the  most  recent  UCSD-processed 
SMEI  data  and  comparisons  with  the  LASCO  CMEs,  the 
HAF  model,  and  ACE  in  situ  data.  In  these  studies  the  SMEI 
images  have  confirmed  that  the  HAFv2  simulations  of  the 
ICME  geometry  are  reasonably  accurate.  However,  in  the 
two  events  of  Howard  et  al/ s  [2007]  paper  and  in  subsequent 
studies,  it  has  been  demonstrated  that  the  HAFv2  model's 
synthetic  skymap  calculation  overestimated  the  ICME 
speed  as  measured  by  SMEI  by  about  a  factor  of  2.  Thus, 
the  SMEI  data  has  already  helped  to  find  and  correct  an 
error  in  the  calculations  in  a  space  weather  model.  (Note 
that  this  error  is  not  present  in  the  shock  arrival  time  fore¬ 
cast  component  of  HAFv2  used  in  the  "Fearless  forecast" 
studies  described  below.) 

[35]  One  goal  is  to  improve  the  operational  solar  wind 
forecasts  at  the  Air  Force  Weather  Agency  (AFWA)  that 
are  based  upon  a  version  of  the  HAFv2  model  which 
became  operational  in  August  2006.  The  technical  objective 
is  to  determine  a  set  of  procedures  that  the  duty  forecaster 
can  use  to  update  or  correct  a  solar  wind  forecast  using 
heliospheric  imager  data  such  as  from  SMEI.  The  following 
scenario  describes  how  this  forecast  adjustment  might  be 
done  in  an  operational  environment. 

[36]  We  assume  that  the  forecaster  has  access  to  a  SMEI 
point-and-click  analysis  tool  and  that  timely  SMEI  images 
are  available  for  analyses.  When  a  solar  event  occurs,  the 
forecaster  will  run  the  HAFv2  model  to  make  an  ICME/ 
shock  prediction.  The  HAFv2  model  also  produces  syn¬ 
thetic  skymaps  of  predicted  ICME  brightness  for  compar¬ 
ison  with  SMEI  images.  The  forecaster  will  use  SMEI 
imagery  to  observe  and  track  the  ICME.  The  forecaster 
will  measure  the  ICME  location  and  speed  using  the  SMEI 
imagery  and  using  the  synthetic  skymaps  produced  by 
HAFv2.  Comparing  the  SMEI-based  and  HAFv2-based 
results,  the  forecaster  will  adjust  a  key  input  to  HAFv2, 
such  as  the  initial  speed  of  the  disturbance  at  the  Sun.  The 
forecaster  will  rerun  HAFv2,  iterating  until  the  observed 
and  forecast  skymaps  match.  The  final  HAFv2  solution 
becomes  the  new  forecast.  When  the  ICME/shock  arrives 
at  (or  does  not  reach)  Earth,  the  forecaster  verifies  the 
forecast  and  updates  the  forecast  skill  statistics. 

[37]  Another  approach  for  comparing  the  forecasting 
capabilities  of  the  HAF  model  and  SMEI  observations  in¬ 
volved  comparing  real-time  "Fearless  Forecasts"  of  shock 
arrival  time  [Fry  et  al.,  2001,  2003]  with  retrospective  SMEI- 
based  forecasts  for  the  20  "halo"  ICME  events  given  by 
Howard  et  al.' s  [2006]  study,  and  computing  the  resultant 
forecast  skill  scores  [see  Fry  et  al.,  2006;  Webb  et  al.,  2006b]. 
An  earlier  study  of  HAFv2  forecasts  of  shock  arrival  times 
(SATs)  during  1997-2002  showed  only  a  moderate  success 
rate  (number  of  correct  forecasts/total  forecasts)  >50%. 
Howard  rt  al.  [2006]  compared  SMEI  ICMEs  associated  with 


Table  1 .  Contingency  Table  for  SMEI  Versus  H  AFv2  Forecasts 


ACE  LI 
Observation 

Forecast 

SMEI- 

Bascd 

Forecast 

HAFv2 

Forecast 

Yes 

No 

Yes 

No 

Yes 

No 

Yes 

hit 

miss 

14 

0 

6 

8 

No 

false  alarm 

correct  null 

6 

0 

0 

0 

Total 

20 

0 

6 

8 

LASCO  halo  CMEs  and  ACE  shocks  at  LI  over  a  1.5-year 
period.  Fry  et  al.  [2006]  determined  the  forecast  hits,  misses, 
false  alarms  and  correct  nulls  for  Howard  et  al.' s  [2006] 
ICMEs  (ex  post  facto  forecasts),  and  compared  these  with 
the  corresponding  real-time  forecasts  by  HAFv2.  For  our 
purpose  a  Hit  =  shock  predicted  within  ±24  hours  of 
observed;  a  Miss  =  shock  observed  but  not  predicted;  a 
False  Alarm  =  shock  predicted  but  not  observed;  and  a 
Correct  Null  =  shock  correctlv  predicted  to  not  reach  Earth. 

[38]  From  these  results  they  formed  a  2  x  2  contingency 
table,  updated  as  shown  in  Table  1.  For  the  14  geomag¬ 
netic  storms  associated  with  the  20  "halo"  ICMEs  given  by 
Howard  et  al.  [2006,  Table  3],  we  issued  8  "Fearless  Fore¬ 
casts"  of  SAT.  Six  of  these  forecasts  were  hits  and  two  were 
misses  (i.e.,  predicted  SAT  minus  observed  SAT  outside  the 
±24-hour  window).  For  6  of  the  storms  where  SMEI 
observed  ICMEs  en  route,  FFs  were  not  issued  because 
the  information  for  the  solar  events  did  not  meet  our  criteria 
for  running  the  model  and  issuing  a  forecast  (as  defined  by 
Fry  et  al.  [2003]).  Had  we  initiated  forecasts  based  upon  the 
SMEI  observations,  6  additional  forecasts  would  have  been 
possible.  These  forecasts  can  also  be  considered  misses, 
bringing  the  total  to  8  misses  and  6  hits.  Thus,  for  this  period 
near-real-time  SMEI  observations,  if  available,  could  have 
nearly  doubled  the  HAFv2  forecast  rate  (i.e.,  the  6  additional 
new  shock  forecasts  versus  the  6  hits  +  2  misses). 

[39]  Models  such  as  the  HAFv2  model  simulations  of  the 
28-29  May  2003  event  can  also  be  used  to  help  explain 
unexpected  features  in  the  SMEI  observations.  In  Sun  etal.'s 
[2008]  simulations  of  this  event,  they  showed  why  this  halo 
ICME  (Figure  1)  appeared  brighter  to  the  east  when  the 
flare  associated  with  the  CME  and  the  assumed  location  of 
the  HAF  energy'  origin  was  west  of  Sun  center.  The  reason 
was  that  a  preexisting  corotating  interaction  region  (CIR) 
likely  contributed  to  the  density  integrated  along  the  line  of 
sight  through  the  eastern  part  of  expanding  ICME  front, 
increasing  the  brightness  of  that  part  of  the  "halo"  (Figure  9). 

[40]  To  summarize  this  section,  we  have  found  that  helio¬ 
spheric  imagers  such  as  SMEI  can  enable  significant 
improvements  in  operational  space  weather  forecasting 
skill.  Specifically  we  compared  forecasts  made  with  the 
HAFv2  model  with  corresponding  SMEI  imagery'.  We 
found  that  SMEI  observations  of  ICMEs,  if  available  in 
near-real  time,  can  be  used  to  make  midcourse  corrections 
to  the  HAFv2  model  forecasts.  We  presented  one  concept  of 
operations  for  adjusting  the  HAFv2  forecast  based  upon 
SMEI  data.  SMEI  observations  and  HAFv2  simulations  are 
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Figure  10.  The  yearly  occurrence  rate  of  the  88  CDAW 
storms.  The  black  curve  is  the  180-day  running  average 
of  daily  sunspot  numbers  in  arbitrary  units.  From 
Zhang  et  al  [2007a,  2007b]. 

complementary  in  that  combining  them  improves  both  the 
accuracy  of  the  HAFv2  model-based  forecasts  and  the 
interpretation  of  the  SMEI  observations. 

4.  Comparing  SMEI  ICMEs  With  the  CDAW 
Intense  Storms 

4.1.  Data  and  Statistical  Analysis 

[41]  Here  we  describe  the  comparison  of  the  SMEI 
ICME  observations  with  the  intense  geomagnetic  storms 
that  were  the  focus  of  the  NASA  LWS  Geostorm  CDAW. 
CDAW  WGl  studied  all  major  geomagnetic  storms  of  cycle 
23  that  occurred  between  January  1996  and  December  2005. 
We  summarize  the  method  of  selection  of  the  storms  as 
discussed  by  Zhang  et  al.  [2007a].  This  10-year  period 
extends  from  the  start  to  late  in  the  declining  phase  of  solar 
cycle  23,  which  had  two  sunspot  maxima  in  2000  and  2001 
(Figure  10).  The  Dst  index  is  a  measure  of  the  strength  of  the 
ring  current  and  is  widely  used  for  measuring  the  intensity 
of  geomagnetic  storms.  A  major  geomagnetic  storm  was 
defined  as  a  minimum  in  the  hourly  Dst  index  falling  below 
Dst  <  -100  nT,  a  commonly  used  threshold  for  major/ 
intense  storms  [e.g.,  Tsurutani  and  Gonzales ,  1997;  Loewe 
and  Pro/ss,  1997].  If  a  period  of  high  activ  ity  showed  multiple 
Dst  <  -100  nT  minima,  CDAW  WGl  arbitrarily  assigned 
these  to  a  single  storm  event  if  the  minima  were  separated 
by  less  than  24  hours,  rather  than  define  each  minimum  as  a 
separate  storm.  CDAW  WGl  identified  88  major  geomag¬ 
netic  storms  during  this  10-year  period,  using  these  selec¬ 
tion  criteria.  These  storms  are  listed  by  Zhang  et  ai  [2007a, 
Table  1],  which  was  corrected  by  Zhang  et  al.  [2007bJ. 

[42]  Since  initiating  operations  in  February  2003,  SMEI 
was  able  to  observe  the  associated  ICMEs  for  most  of  the 
CDAW  intense  storms.  Here  we  describe  the  SMEI  obser¬ 
vations  and  analyses  for  the  26  CDAW  storms,  CDAW 
numbers  63  -  88,  that  occurred  during  the  SMEI  observa¬ 
tions  from  May  2003  to  October  2005  In  addition  for 
completeness,  we  include  the  two  other  intense  storms  that 
occurred  after  that  period  through  the  end  of  2007  for  a  total 
of  28  storms.  The  solar  and  IP  drivers  of  these  latter  storms 
were  unambiguous. 


[43]  These  28  storms  are  listed  in  Table  2,  along  with 
details  of  the  SMEI  observations  and  of  the  CDAW-iden- 
tified  solar  and  interplanetary  (at  1  AU)  sources  of  the 
storms.  Details  of  the  search  and  identification  process 
are  given  by  Zhang  et  al.  [2007a]  and  will  not  be  repeated 
here.  For  our  study  we  assume  that  these  identifications 
were  correctly  made  and  we  will  use  the  timing  and  location 
information  as  given  by  Zhang  et  al.  [2007a,  Table  1].  We 
merely  note  that  the  identification  procedure  used  by 
CDAW  WG  1  focused  on  identifying  the  chain  of  activity 
preceding  a  storm  starting  with  the  solar  CME  and  associ¬ 
ated  surface  activity  and  then  the  interplanetary  counter¬ 
parts  of  this  ICME  at  1  AU,  especially  considering  evidence 
of  shocks,  magnetic  clouds  and  other  ICME  material.  The 
primary  data  sets  used  were  from  the  SOHO  LASCO 
coronagraphs  and  EIT  instrument  at  the  Sun,  and  from 
the  ACE  spacecraft  at  LI  and  the  Wind  spacecraft  near 
Earth.  In  addition,  we  note  that  the  SMEI  observations 
discussed  herein  were  used  by  the  working  group  to  help 
track  ICMEs  to  larger  distances  from  the  Sun  than  is 
possible  with  LASCO  and,  therefore,  to  aid  in  the  identifi¬ 
cation  and  timing  of  the  Earth  arrival  of  the  shock  and/or 
ICME,  and  the  storm  onset. 

[44]  Table  2  presents  the  basic  data  and  some  results  for 
our  study  of  SMEI's  observations  of  ICMEs  associated  with 
the  CDAW  storms.  The  first  three  columns  of  Table  2  give 
for  each  storm  its  CDAW  number,  the  peak  intensity  in  - 
Dst,  and  its  estimated  onset  (first  row)  and  peak  (next  row) 
times.  The  Dst  data  were  obtained  using  the  Kyoto  geo¬ 
magnetic  data  at  http://swdcdb.kugi.kyoto-u.ac.jp/dstdir/ 
index.html.  The  storm  values  through  2003  were  based  on 
the  final  Dst  index,  whereas  those  in  2004  and  2005  were 
based  on  the  provisional  Dst  index,  so  that  the  final  values 
may  change  slightly.  We  estimated  onset  times  from  the 
hourly  index  based  on  the  hour  when  Dst  began  a  system¬ 
atic  drop  starting  at  ^-40  nT  until  it  reached  its  minimum 
intensity.  The  next  eight  columns  describe  the  SMEI  data 
for  each  ICME  that  we  associated  with  each  storm.  The  last 
nine  columns  are  from  the  CDAW  Table  1  and  describe  the 
properties  of  the  solar  source(s)  and  the  properties  of  the  IP 
source(s)  that  the  CDAW  working  group  decided  directly 
"drove"  the  geomagnetic  activity. 

[45]  For  the  SMEI  data,  in  Table  2,  in  the  Time  First 
Observed  column  two  dates  (and  day  of  year)  and  times 
are  given  when  SMEI  first  observed  the  ICME  (top)  and 
the  estimated  time  (hour)  when  it  left  the  Sun  (next  row). 
This  "onset"  time  uses  a  linear  extrapolation  of  data 
points  on  an  5  versus  time  plot  back  to  0  elongation, 
i.e.,  the  Sun.  The  Storm/CME  Association  column 
indicates  whether  or  not  we  considered  the  SMEI 
transient  to  be  associated  with  the  storm  (top)  and/or 
with  the  LASCO  CME(s)  (next  row)  in  the  Solar  col¬ 
umn.  In  the  Aurora  Date  column  we  give  the  date  and 
time  when  SMEI  first  observed  auroral  light  associated 
with  the  storm.  In  Table  2  the  f^iff  Au-FO  column  gives 
the  time  difference  in  hours  between  the  ICME  first 
observed  ("FO")  time  and  the  auroral  ("Au")  onset 


13  of  24 


S05002 


YVEBB  ET  AL.:  GEOEFFECTIVE  ICMES  BETWEEN  SUN  AND  EARTH 


S05002 


s 

p 

X 

a> 

x 


0> 

2 


<N 

a; 

2 

rj 

H 

2 

at 

£ 

n 

c/3 

Ui 

S 

u 

2 

£ 

LD 

73 

C 

a 


o 

2 

£ 

< 

a 

u 

2 

3- 

E 

ra 

5/5 

rsi 

2 

2 

f- 


c 

o  £  * 

*X3  C 
"J  c 

g  *  Sfi 


a.  m 

c/3  aT 
73 


Or? 

-?<  V)  £ 


3  £ 
£<  3 
"3  X  C 
**(^X 


*  c 

73  D  C  6  « 

*  <  3  H 

>  73  ! 


-  .g  tj  i 

I s 


E  e 

CU  -J3 

ilE 

i:  3 ! 

uj  x 
20 
cx 


a 

a  g 
OE^ 
*9  H  fc 

&T.5 

i3  C 
<  *" 


2  £ 
y  a 
E’g 
S  I 

C/3  < 


-n  X  x 

a  *  2 

rj  O  £  p1 

J 

“c5 


E  .* 

H  3  _ 

73  H 1 
>  a  4J  ^ 

1  „  55 
i  £  c 
a  C 
C 


c  C 


2  E 

<  1 
C  5 

U  </> 


S  § 


o  , 

(N 
V  ' 


2 

rJ 

£ 

rj 

£ 


a 

<N 

>; 

P3 

2 


<N 


bC 

3 

< 


5 

0 

O' 

<N 

? 


s  s 


<b 

2 


2 
r  • 
2 


3  CL 


r. 

u 

r. 

*  >. 

§ 

H 


3 

< 


> 

C 

2 


t: 

o 


££>-  >• 


O'  sC 
<^3  (7 

O'  ad 

?— i  r~< 

>  > 
c  c 
2  2 
OH  o 

ri  CM 


I  § 

v  D- 

!  & 
z  2 
So 

y  <N 


>  > 


£  S*s 


r- 

'-'W  <!  C 

^si 

ac 

rsi 


a 


C  i 


« 9 8 Zi § 5  2  =  3 

■  Sf5§t23t!8jf  o  ^  = 

>  *-s.t<0  C’O  oSv'i 

S* 


£  2  2, 


rT  2  pj  oe.  3  wy  nr. 
^  IT  <  r-  CM 


as  2 

3  3^ 
<  <0 
t>.  x  2s, 


s 


If 

O  t> 


<N  £4 

?  ct 


'tOOMTHNn 
CO  r-  cr>  r*  m  in 
vC  ^  r*'  N  (N 
r-  O  (S  O  H  rs] 


Cl.  Cl.  Cl,  CL  3  3 

«<<M 
N  r  r  (N  N 


a 


a 

< 


er.  > 

a  £  o 

^  3  X 

£  . .5  «* 

5  B  E 
<  o  o 


§s>.- 
'  u  a  > 


>•  >• 


g 

g 

8 

II 

gg 

8  8 

gg 

gg 

§8 

O 

- 

S 

r4 

CM  rr 

t~<  T— 1 

ur.  x 

r~<  *—> 

X  r- 

o 

as 

i—  <N 
O  r- 

tf 

c 

u 

c 

tf 

G 

>  > 
c  c 

2  2 

c  c 

n 

-a  -a 
ti  u. 

Cl.  Cl. 

<  < 

3  3 

3  3 

o 

o 

rr> 

o 

rri 

o  c 

CM  <N 

aa 

r~>  1—* 

r~i  t— < 

X.  -1* 

aa 

aa 

<r, 

O' 

<N 

X 

X 

rr, 

Tf 

r~i 

O 

PH 

c 

Tf 

5 

3" 

X 

O 

r>. 

N. 

x. 

N. 

*t 

K. 

£  JQ 

£5 


"Si 

8-3 

c 

c 


V) 


cr 


ai  — 
ai 

w  C 

a  ai 


a;  v, 
a»  £ 

r*  *9 
■3 

£  E 
£  S 

C-  X 
£ 

-a  o 
3  X 

iA 

«  £ 

2  x 

c 

73  v3 

X.E 

3  « 
o  E 

C  -3 
a»  73 

Jr  C 


o 

>> 

rj 

73 


3  9S 

_,  c 

73  c 
a» 

N  73 
'C  aj 

2  * 
5  E 

5  X 


a  . 

£  g 


r-  ^  TJ 

3  >  = 

2  >  x 


14  of  24 


S05002 


WEBB  ET  AL.:  GEOEFFECTTVE  ICMES  BETWEEN  SUN  AND  EARTH 


S05002 


(top)  and  storm  ("St")  onset  (next  row)  times.  Finally, 
we  give  the  estimated  span  of  the  ICME  and  its 
estimated  elongation  at  the  time  of  storm  onset,  both 
in  degrees.  These  data  and  results  will  be  discussed 
later.  The  Predicted  3-D  and  tdiff  St-Au  columns  are 
also  discussed  later. 

[46]  In  Table  2,  in  the  Solar  column  is  listed  the  solar 
source  type  as  S,  single  CME  (18  events);  M,  multiple 
CMEs  (6);  and  CH,  coronal  hole  (4).  Considering  the 
properties  of  solar  sources,  the  Onset  Date  and  Time 
column  gives  the  first  appearance  of  the  CME  in  the 
LASCO  C2  coronagraph,  the  CME  Speed  column  gives 
the  average  speed  of  the  CME  through  a  linear  fit  in  the 
LASCO  C2/C3  fields  of  view,  and  the  Width  column  gives 
the  apparent  angular  span  of  the  CME  in  the  plane  of  the 
sky  measured  in  the  C2  field  of  viewT.  These  values  were 
generally  obtained  from  the  online  LASCO  CME  catalog  at 
http://cdaw.gsfc.nasa.gov/CME_Iist/  [Yashiro  et  al,  2004]. 
Source  Location  column  gives  the  heliographic  coordi¬ 
nates  of  the  associated  surface  activity.  This  generally 
corresponds  to  the  Ho  flare  location  reported  by  the 
NOAA  Space  Weather  Prediction  Center  (SWPC).  When 
no  Ha  flare  location  was  reported,  SOHO  EFT  images 
were  used  to  measure  the  source  coordinates,  given  by 
the  location  of  any  brightening  or  the  dimming  regions 
observed.  The  surface  identification  of  the  candidate  CME 
for  storm  76  was  considered  inconclusive,  however  we  did 
not  consider  that  storm  because  of  poor  SMEI  data. 
Because  of  SOHO  data  gaps,  the  solar  source  for  storm 
72  was  identified  with  a  major  long-duration  solar  flare  in 
the  GOES  SXI  imager  and  a  halo  CME  observed  a  few 
hours  later  in  LASCO  C3. 

[47]  In  Table  2,  for  the  6  M-type  events,  numbers  63,  64, 
77,  78,  79,  86,  there  are  multiple  rowrs  for  each  event  listing 
the  multiple  CMEs  that  may  have  contributed  to  the 
observed  1  AU  solar  wind  structures.  In  each  case,  the  first 
row  indicates  what  is  considered  the  "principal"  solar 
driver.  In  the  case  of  the  4  C-type  events,  the  definitions 
of  the  parameters  in  the  solar  source  columns  are  slightly 
different  because  of  the  different  nature  of  the  source.  The 
time  in  the  Onset  column  is  the  central  meridian  transit 
time  of  the  centroid  of  the  associated  coronal  hole  mea¬ 
sured  from  EIT  images.  The  time  is  followed  by  "(CH)"  in 
order  to  emphasize  that  this  does  not  refer  to  a  CME  source. 
The  data  for  the  storms  driven  by  CH-CIRs,  numbers  65,  71, 
81,  87,  are  italicized;  these  events  were  not  used  in  this 
study. 

[48]  For  the  properties  of  the  IP  sources,  in  Table  2,  the 
IP  column  summarizes  the  solar  wind  components  at  1  AU, 
in  time  order,  that  CDAW  WG1  identified  as  driving  the 
storm.  The  primary  drivers  were  considered  to  be  periods 
of  southward  magnetic  field  either  in  shock-related  mate¬ 
rial  or  in  the  ICME.  The  structures  in  IP  are  SH,  shock 
sheath;  ICME,  interplanetary  CME;  MC,  magnetic  cloud. 
(M)  indicates  multiple  structures  of  that  type.  A  hyphen 
indicates  an  interaction  between  two  structures,  in  partic¬ 
ular,  preceding  ICME  (PICME)-SH  (numbers  64  and  78) 


denotes  a  shock  propagating  through  a  preceding  ICME. 
The  Shock  and  ICME  Onset  columns  show  the  date  and 
time  of  the  ICME-driven  shock  from  ACE  or  Wind  and  the 
onset  time  of  the  ICME,  respectively.  For  number  64  the 
first  row  is  an  estimate  for  the  shock  arrival  time.  The  ff 
Au-Sh  column  will  be  discussed  later. 

[49]  Next  we  discuss  several  examples  illustrating  the 
SMEI  geoeffective  ICME  observations  and  how  they  related 
to  solar  wind  conditions  in  front  of  Earth  and  the  storm 
effects  at  Earth.  The  first  example  was  introduced  earlier. 
This  was  the  first  Earth -directed  ICME  detected  and  ana¬ 
lyzed  in  the  SMEI  data,  occurring  on  28-29  May  2003, 
5  months  after  observations  began  [Tappin  et  al ,  2004; 
Jackson  et  al. ,  2007]  (Figures  1  -3).  This  fast  (1000  km/s)  event 
was  detected  at  —30°  e,  about  15  hours  before  it  passed 
over  the  Earth  on  29-30  May  2003,  causing  a  major 
geomagnetic  storm.  The  ICME  erupted  from  the  Sun 
following  two  bright  (X-class)  flares  from  an  active 
region  near  Sun  center.  These  yielded  two  nearly  si¬ 
multaneous  halo  CMEs  in  LASCO  which  SMEI  later 
sawf  as  three  contiguous  arcs  together  covering  over 
150°  of  sky  (Figure  lb). 

[50]  The  elongation-time  plot  of  the  leading  edges  of  the 
LASCO  CME  and  SMEI  ICME  features  are  shown  in 
Figure  2  (top),  along  with  time  markers  at  90°  £  indicative 
of  the  ATs  at  1  AU  (from  Table  2)  of  solar  wind  shocks 
(asterisks),  the  SMEI  aurora  onset  time  (crosses),  and 
the  storm  onset  and  peak  times  (the  two  bars,  respec¬ 
tively).  The  bottom  plot  is  in  "distance"  versus  time 
wherein  we  converted  the  LASCO  CME  and  SMEI 
ICME  £-time  data  by  assuming  that  the  CME  was 
launched  from  the  CDAW-identified  source  location 
and  propagated  radially  from  the  Sun.  This  method  is 
discussed  later  in  section  4.2.1.  Here  the  time  markers 
are  accurately  located  at  1  AU.  These  data  were  then 
used  to  calculate  the  arrival  times  of  the  ICME  struc¬ 
tures  at  1  AU  based  on  linear  fit  extrapolations  of  the  3 
separately  tracked  SMEI  features  (the  two  later  ATs  are 
given  in  the  Predicted  3-D  column  of  Table  2).  This 
shows  that  two  of  the  SMEI  features  might  have  been 
associated  with  the  two  shocks,  with  the  SMEI  aurora 
onset  time  lying  between  them.  The  third,  and  best 
observed  SMEI  arc  was  closely  associated  with  the 
storm  onset  itself.  This  timing  is  confirmed  in  the 
solar  wind  and  geomagnetic  data  plot  of  Figure  3.  Key 
solar  wind  plasma,  composition  and  IMF  parameters 
are  shown  in  these  plots  (see  also  Figures  5f-5p  and 
12)  as  well  as  the  Dst  and  Kp  geoindices.  The  vertical 
red  line  denotes  the  SMEI  aurora  onset  time  between 
the  two  shocks,  and  about  7  hours  before  the  Dst 
storm  onset  (however,  note  that  the  3-hour  Kp  index 
is  enhanced  earlier).  CDAW  WGl  determined  that 
this  storm  was  caused  by  southward  IMF  mostly 
between  the  two  shocks,  with  multiple  and  com¬ 
pressed  ICME  material  involved.  There  was  clear 
ICME  material  and  strong  IMF  following  the  second 
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shock  but  the  IMF  was  strongly  northward,  minimizing 
any  coupling. 

[si]  During  two  weeks  in  late  October  and  early  No¬ 
vember  2003,  a  series  of  large  solar  events  led  to  high  levels 
of  energetic  particles  in  geospace  and  produced  overlap¬ 
ping  large  geomagnetic  storms  on  28-30  October.  Erup¬ 
tions  came  from  three  active  sunspot  regions.  The  first  rotated 
into  view  on  19  October  and  over  the  next  two  weeks  the 
three  regions  unleashed  11  major,  X  class  flares,  an  unpre¬ 
cedented  production  rate  for  such  powerful  events.  The 
three  most  geoeffective  flares  and  halo  CMEs  occurred  on 
22, 28,  and  29  October.  The  weaker  22  October  ICME  caused 
only  a  brief  geomagnetic  storm  that  was  terminated  by 
strong  northward  fields  within  the  ICME.  However,  the 
28  and  29  October  flares  had  peak  fluxes  as  measured  by 
the  NOAA  GOES  satellites  (X17  and  X10)  that  were  among 
the  largest  ever  recorded,  and  were  associated  with  large 
particle  events  at  Earth  peaking  early  and  late,  respectively, 
on  29  October.  The  shock  waves  from  these  events  arrived 
at  Earth  only  19  hours  after  leaving  the  Sun,  among  the 
fastest  events  ever  observed  since  the  dawn  of  the  satellite 
age  [Cliver  et  al,  1990].  These  shocks  and  their  trailing 
ICMEs  drove  two  intense  storms,  CDAW  storms  67  and 
68,  with  peak  Dst  values  both  on  30  October  of  353  and 
383  nT,  respectively. 

[52]  SMEI  observed  the  22  and  28  October  halo  ICME 
events,  again  starting  at  distances  about  1/3  of  the  way  from 
the  Sun  to  Earth  and  21  and  10  hours,  respectively,  before 
Earth  arrival.  However,  the  second  energetic  event,  on  29- 
30  October,  was  one  of  the  7  CDAW  storms  for  which  SMEI 
had  no  or  insufficient  data.  This  was  because  the  intense 
particle  bombardment  and  aurora  from  the  first  SEP  event 
and  storm  knocked  out  the  instrument  through  most  of 
30  October  and  into  31  October!  For  the  28-29  October 
storm  (67,  plot  not  shown),  SMEI  observed  two  ICME  arcs; 
the  d-t  plot  yielded  ICME  ATs  at  1  AU  between  the  shock 
AT  and  storm  onset  and  after  storm  onset,  respectively 
(Appendix  A). 

[53]  Figures  5  and  6  show  SMEI  and  solar  wind  data  during 
the  period  in  November  2004  involving  the  two  CDAW 
storms  77  and  78.  The  solar  wind  plot  in  Figure  5  shows 
that,  although  there  were  a  number  of  shocks  during  this  4-5 
day  period,  there  were  clearly  two  long  ICME  intervals 
(Figure  5o)  associated  with  the  strong,  southward  IMF 
driving  the  storms.  These  intervals  are  denoted  by  magnetic 
clouds  (Figures  5e-5g),  low  Tp  and  density  (Figures  5h  and 
5i),  and  high  He  and  07/06  content  (Figures  5k  and  51).  The 
storms  were  driven  (1)  in  CDAW  77  by  high  speed  and 
southward  fields  associated  with  multiple  shocks  and  the 
first  cloud  and  (2)  in  CDAW  78  a  shock  propagating  through 
this  preceding  ICME  (note  the  Dst  dip  labeled  'D')  followed 
by  another  shock  and  the  second  cloud.  SMEI  observed  at 
least  two  large  arcs  or  partial  halos  east  and  southeast  of 
the  Sun  from  5  to  7  November  that  may  have  merged  near 
Earth  (Figures  5a -5c).  Two  fainter  arcs  were  seen  off  the 
west  limb  on  8-9  November  (Figure  5d)  that  were  possibly 
associated  with  dip  D  and/or  storm  78.  The  vertical  red 


CDAW  STORM  NUMBER  73 


Figure  11.  Example  of  an  Earthward  partial  halo 
ICME/ICME  and  CDAW  storm  73.  (a)  LASCO  C2  differ¬ 
ence  image  on  20  July  2004  at  1454  UT.  (b)  SMEI  differ¬ 
ence  images  on  22  July  at  0530  and  1040  UT  showing  the 
halo  structures  (arrows).  The  storm  onset  was  on  22  July 
at  —1500  UT  (Figure  8).  (c)  Distance-time  plot  and 
(d)  elongation -time  plot.  Symbols  are  as  in  Figure  2. 


bars  on  Figures  5e-5o  indicate  the  times  when  the  aurora 
associated  with  the  two  storms  were  first  observed  by 
SMEI.  The  ICME  arcs  leading  edges  were  at  105°  and  85° 
£,  respectively,  at  these  times  that  were  each  about  half 
a  day  before  the  storms'  main  phases.  This  shows  the 
lead  times  that  a  SMEI-type  imager  could  provide. 

[54]  Figure  6  shows  the  e-f  and  d-t  plots  for  this  period, 
with  the  LASCO  CME  and  SMEI  ICME  fronts  and  the 
time  markers  at  1  AU.  Both  plots  show  the  1  AU  AT  of  the 
first  ICME  between  the  first  shock  AT  and  the  SMEI 
aurora  onset  time,  and  the  second  ICME  AT  coincident 
with  another  shock  and  several  hours  before  the  Dst  dip, 
D,  that  preceded  storm  78.  Thus,  the  second  SMEI  ICME 
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Figure  12.  ACE  solar  wind  and  storm  data  for  22-23  July  2004  ICME  and  CDAVV  storm  73.  Plot  is 
as  in  Figure  3.  The  vertical  red  bar  indicates  the  onset  time  of  bright  aurora  in  SMEI  (arrow)  when 
the  ICME  front  was  at  an  elongation  of  60°.  The  vertical  green  line  indicates  the  shock  AT  at  ACE 
(see  text).  ACE  plot  courtesy  of  I.  Richardson. 


may  have  been  associated  with  this  dip.  Although  of  short 
duration,  this  dip  was  still  significant  (peak  =  223  nT)  and 
likely  driven  by  the  compression  due  to  the  second  shock 
and  ICME  running  into  the  first. 

[55]  Our  final  example  involves  CDAW  storm  73  in  July 
2004.  A  large-scale  eruption  occurred  over  the  northwest 
part  of  the  Sun  on  20  July  leading  to  a  partial  halo  CME  in 
LASCO  over  the  entire  northern  hemisphere  (Figure  11a). 
Later,  on  21-22  July  SMEI  observed  the  ICME  as  several 
expanding  arcs  covering  the  northern  ecliptic  hemisphere 
reaching  spans  of  120°  (Figure  11c).  Thed-f  plot  of  Figure  11c 
shows  how  well  the  LASCO  and  SMEI  points  line  up  for  this 
event,  and  that  the  SMEI  predicted  AT  at  1  AU  is  nearly 
simultaneous  with  the  SMEI  aurora  onset  time  at  22  July  at 


1500  UT.  This  time  is  also  about  7  hours  between  the  shock 
AT  and  storm  onset  (see  the  solar  wind  data  for  this  event  in 
Figure  12),  illustrating  the  general  pattern  found  for  all  the 
CDAW  storms  (see  below).  CDAW  WGl  concluded  that 
the  storm  was  driven  by  southward  field  and  high  speeds  in 
the  ICME  following  the  shock  and  its  sheath  region. 

[5f>]  Next  we  provide  some  statistical  results  derived  from 
the  tabular  data.  Figure  10  presents  the  annual  occurrence 
rate  of  the  88  intense  CDAVV  storms  with  sunspot  numbers 
superposed  on  it.  We  see  that  the  storm  rate  peaks  with  the 
sunspot  numbers  in  2000-2001,  but  that  there  are  many 
storms  during  the  rise  and  declining  phases  of  the  solar 
cycle.  The  storms  of  this  study  include  those  in  the  last 
3  bins,  2003-2005  plus  the  two  additional  events  through 
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2007,  all  occurring  during  the  declining  phase  of  this  cycle. 
Of  the  28  intense  storms  during  this  period,  7  have  either  no 
or  inadequate  SME1  data  coverage  preceding  the  storm  peak, 
leaving  21  events  for  further  study.  Three  of  these  21  storms 
were  identified  during  the  CDAW  as  being  driven  solely  by 
CIRs  associated  with  solar  coronal  holes  (numbers  65,  71,  87). 
The  CIR-driven  CDAW  storms  were  discussed  by  Richardson 
etal.  [2006].  Since  it  has  not  yet  been  demonstrated  that  SMEI 
can  detect  CIRs,  we  will  not  consider  these  storms. 

[57]  This  leaves  18  remaining  storms  with  adequate  SMEI 
coverage  and  for  which  CDAW  WG1  identified  the  solar 
and  IP  sources.  Of  these  18  storms,  12  were  preceded  by  a 
definitely  associated  SMEI  ICME,  4  by  a  likely  associated 
ICME  and  2  had  no  associated  ICME.  Thus,  we  can  con¬ 
clude  that  for  89%  (16/18)  of  the  intense  storms  during  this 
period,  SMEI  likely  detected  and  tracked  the  lCMEs  that 
caused  the  storms. 

[58]  The  storm-SMEI  ICME  associations  were  done  in¬ 
dependently  of  knowledge  of  the  solar  and  IP  sources.  Thus, 
we  can  evaluate  whether  or  how  well  the  16  SMEI  ICMEs 
associated  with  the  storms  could  be  identified  with  the 
CDAW-source  CMEs.  In  these  16  cases,  the  SMEI  ICMEs 
and  LASCO  CMEs  were  definitely  associated  in  9  cases, 
likely  associated  in  5,  and  questionable  in  2.  Thus,  in  87.5% 
(14/16)  of  the  storm-SMEI  ICME  events,  SMEI  likely 
detected  the  geoeffective  ICME  observed  by  LASCO  leav¬ 
ing  the  Sun.  On  the  other  hand,  since  all  the  storms  had 
reasonably  reliable  solar  CME  identifications,  this  suggests 
that  SMEI  did  not  clearly  detect  the  appropriate  ICME  in 
4  cases,  the  two  questionable  cases  (numbers  85  and  86)  and 
the  two  other  cases  with  no  SMEI  ICME  (numbers  64  and 
74).  However,  in  event  85  CDAW  WGl  concluded  that  the 
solar  source  was  inconclusive,  and  there  were  either  SOHO 
or  SMEI  data  gaps  during  the  other  events.  Thus,  even  in 
some  or  all  of  these  4  "no  detection"  cases,  SMEI  might 
have  detected  the  geoeffective  ICME  given  full  coverage. 

[59]  As  discussed  earlier,  one  of  the  unexpected  discov¬ 
eries  of  SMEI  is  that  it  is  able  to  detect  auroral  light  above 
its  orbital  height  during  geoactivity  periods  when  Kp  >  ~4 
[Mizuno  et  al,  2005].  The  intense  storms  considered  here 
typically  reach  higher  levels  of  Kp,  and  the  Kp  =  4  level 
corresponds  to  the  Dst  =  -40  nT  level  we  used  to  define 
the  onset  of  each  storm  [Loezve  and  Prolss,  1997].  Thus,  the 
SMEI  sky  maps  could  readily  be  used  to  determine  the 
onset  time  of  each  of  the  CDAW  storms,  even  in  those  cases 
where  SMEI  did  not  detect  the  oncoming  ICME,  or  there 
was  a  CIR  source.  In  the  td iff  St-AU  column  of  Table  2  we 
show  for  each  storm  the  time  difference  in  hours  betw  een 
when  the  auroral  light  associated  with  the  storm  first 
appeared  on  the  SMEI  sky  maps  and  the  storm  onset  time 
defined  above.  Excluding  lower  limits,  the  mean  of  the 
algebraic  differences  is  +4.1  hours  (10  events;  range  = 

2.5  to  +8.15  hours).  Historically  storm  sudden  commence¬ 
ments  (SSCs)  tend  to  precede  the  main  phase  of  storms. 
Since  SSCs  are  associated  with  shock  waves  in  the  solar 
wind,  it  is  of  interest  to  compare  the  timing  between  the 
CDAW  shock  times  observed  by  either  ACE  or  Wind  and 


the  onset  time  of  the  SMEI  aurora.  So,  in  the  Shock  column 
we  list  the  time  difference  in  hours  between  when  SMEI 
first  observed  the  aurora  and  the  1  AU  shock  time.  Exclud¬ 
ing  lower  limits,  the  mean  of  the  algebraic  differences  is 
+3.9  hours  (14  events;  range  =  -22.9  to  +15.4  hours).  Thus, 
on  average  the  shock  arrived  ~4  hours  before  the  SMEI 
auroral  onset,  which  in  turn  preceded  the  storm  onset  by 
^4  hours.  These  time  differences  probably  have  uncertain¬ 
ties  of  1-2  hours,  because  SMEI's  time  resolution  is  dic¬ 
tated  by  its  102-min  orbit  and  Dst  is  a  1-hour  averaged 
index.  We  have  also  not  corrected  for  the  time  a  shock  takes 
to  travel  the  distance  between  the  LI  point  where  ACE  is 
and  Earth.  For  typical  solar  wind  speeds  this  travel  time  is 
approximately  40  min. 

h>o]  Thus,  the  typical  pattern  for  these  intense  events 
was  (1)  arrival  of  the  IP  shock  followed  ^4  hours  later  by 

(2)  the  SMEI  auroral  onset  followed  --4  hours  later  by 

(3)  the  storm  onset  ( Dst  <  -40  nT  criterion).  We  can  think 
of  two  reasons  why  the  SMEI  aurora  precedes  storm  onset 
by  a  few  hours.  First,  it  is  known  that  IP  shock-driven 
SSCs  usually  are  associated  with  a  sudden  brightening  of 
the  aurora  [ Craven  etal.,  1985;  Egeland  et  al,  1994;  W.  Burke, 
private  communication,  2008].  Second,  this  suggests  that 
SMEI  may  be  sensitive  to  even  lower  auroral  light  levels 
than  the  Dst  <  -40  nT  level,  and,  therefore,  might  begin  to 
detect  the  aurora  before  the  -40  nT  level  is  reached.  Of 
course,  the  Dst  index  is  defined  from  low-latitude  geo¬ 
magnetic  observatories  that  sample  effects  in  the  ring 
current  at  altitudes  of  several  hundred  km.  Since  SMEI 
views  light  emission  above  heights  of  a  few  thousand  km, 
it  is  seems  remarkable  that  the  correspondence  with  the 
Dst  index  is  so  good. 

[6i]  In  Table  2,  in  the  fdiff  Au-FO  column  (top  row)  we 
show  the  time  difference  between  when  SMEI  first  ob¬ 
served  the  ICME  and  when  the  SMEI  aurora  first 
appeared.  The  next  row  of  the  Au-FO  column  gives 
the  time  difference  between  the  SMEI  ICME  first  observed 
time  and  the  storm  onset  time.  The  means  of  these  two 
time  differences  are  23.2  hours  (11  events;  range  =  7- 
36  hours)  and  27.5  hours  (15  events;  range  =  14-44  hours), 
respectively.  (The  difference  in  these  two  mean  values,  of 
course,  is  the  4-hour  average  difference  discussed  above 
between  the  SMEI  aurora  and  storm  onset  times.)  From  a 
space  weather  standpoint,  the  most  important  conclusion 
from  these  timing  studies  is  that  SMEI  could  detect  and 
begin  to  track  the  geoeffective  ICME  an  average  of  28 
hours  before  it  "struck"  Earth,  in  terms  of  the  Dsf- defined 
storm  onset.  Thus,  this  represents  an  average  lead  time  or 
warning  time  for  a  storm  which  could  be  provided  by  a 
similarly  designed  heliospheric  imager  with  real-time  data 
latency  capability.  Since  the  average  transit  time  of  the 
ICMEs  (storm  onset  -  CME  onset)  was  54  hours,  this 
suggests  that  SMEI  could  begin  to  track  them  halfway 
between  the  Sun  and  Earth.  The  range  in  this  warning 
time  was  0.6- 1.8  days,  depending  on  such  factors  as 
ICME  speed,  brightness,  data  coverage,  etc.  (the  range 
of  ICME  transit  times  was  1.15-3.2  days). 
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[62]  In  Table  2 ,  the  Span  column  shows  the  estimated 
maximum  angular  span  in  degrees  of  the  geoeffective 
ICME  structure  observed  by  SMEI.  The  angular  span  is 
the  difference  in  Position  Angle  between  the  observed 
sides  of  the  ICME  structure.  This  is  often  a  lower  limit 
because  of  interference  by  particles,  auroral  light,  and 
saturated  or  shutter-closed  frames.  The  mean  span  of 
78.6°  (11  events;  range  =  30-130°)  is  typical,  or  slightly 
larger  than  those  of  the  Earthward,  arc-like  events  dis¬ 
cussed  in  section  2.2.  This  large  size  is  not  unexpected  for 
a  plasma  cloud  that  is  observed  as  it  is  traveling  near  the 
Earth. 

[63]  In  Table  2,  the  Elongation  column  gives  the  esti¬ 
mated  elongation  of  the  leading  edge  of  the  ICME  ob¬ 
served  by  SMEI  at  the  time  of  storm  onset.  The  mean  £  = 
78.7°  (15  events;  range  =  45-105°).  Calculations  suggest 
that  such  large  elongation  values  are  expected  for  solar 
transients  traveling  to  1  AU  and  beyond  and  to  within 
the  Earth's  vicinity  [e.g.,  Kahler  and  Webb,  2007].  How¬ 
ever,  for  events  directed  Earthward,  the  observations  are 
complicated  by  the  varying  geometry  and  size  of  the 
ICME  as  well  as  by  its  brightness  which  is  affected  by 
Thomson  scattering  [Vourlidas  and  Howard ,  2006].  That 
the  geoeffective  ICMEs  have  both  large  spans  and  large 
elongations  at  storm  onset  does  not  necessarily  tell  us 
much  about  their  inherent  size  or  geometry,  only  that 
they  are  being  viewed  nearby,  i.e.,  when  close  to  Earth. 
Their  arc-like  nature  suggests  that  we  are  viewing  the 
emission  across  a  spherical  surface  but  we  do  not  know 
whether  the  true  "nose,"  or  radial  center  of  curvature,  is 
along  that  surface.  See  the  appendix  of  Kahler  and  Webb 
[2007]  for  a  discussion  of  geometry  versus  distance 
measurements  of  ICMEs. 

4.2.  Forecasting  Implications  of  the  SMEI 
Results 

4.2.1.  Using  SMEI  Data  to  Predict  Geoeffective 
ICME  Arrival  Times 

[64]  We  can  now  use  the  results  of  the  SMEI-CDAW 
study  to  update  and  extend  Howard  et  al.'s  [2006]  result  dis¬ 
cussed  in  section  3.2,  specifically  to  produce  a  new  version 
of  the  time  difference  results  that  we  produced  from  their 
data  (Figure  8).  Their  study  compared  SMEI  ICMEs  asso¬ 
ciated  with  LASCO  halo  CMEs  and  ACE  shocks  at  Ll  over  a 
1.5-year  period.  Our  study  began  with  the  CDAW  storms 
and  their  identified  CMEs  and  source  regions  and  searched 
for  associated  SMEI  ICMEs  from  2003  to  2005  (2007).  De¬ 
spite  these  differences  there  is  an  overlap  of  8  events 
between  the  two  studies  (CDAW  63,  64,  66,  67,  69,  70,  72, 
and  73).  Howard  et  al  [2006]  computed  distance- time  pro¬ 
files  using  the  point  P  approximation  to  fit  the  SMEI  ICME 
data  and  predict  their  times  of  arrival  at  Earth.  They  then 
compared  these  predicted  ICME  arrival  times  with  the 
actual  shock  times  at  Ll,  and  we  resorted  their  results  to 
produce  the  time  difference  histogram  of  Figure  8. 

[65]  In  this  study  as  stated  earlier,  we  converted  the 
LASCO  CME  and  SMEI  ICME  -time  data  to  distance¬ 


time  by  assuming  that  the  ICME  was  launched  from  the 
CDAW-identified  source  location  and  then  propagated 
radially  away  from  the  Sun.  This  is  related  to  the  "fixed 
4>"  method  described  by  Howard  et  al.  [2007]  and  Kahler 
and  Webb  [2007].  Figures  2  (bottom),  6  (bottom),  and  11c 
are  examples  of  these  distance-time  plots.  We  made  no 
assumptions  about  the  speed  or  acceleration  of  the 
ICME.  We  do  implicitly  make  the  assumption  that  both 
LASCO  and  SMEI  are  tracking  the  same  front  of  mate¬ 
rial  as  it  moves  out.  Despite  the  difficulties  that  can  be 
raised  with  each  of  these  assumptions,  our  goal  here  is 
merely  to  estimate  the  uncertainties  in  arrival  times  that 
such  straightforward  assumptions  would  yield. 

[66]  To  determine  the  distances  it  was  necessary  to 
attempt  to  correct  for  the  projection  effects  that  are  typical 
of  solar  and  heliospheric  images  and,  thus,  to  convert  the 
angular  distance  to  a  "true"  3-D  distance.  This  was  done 
using  the  following  equation  from  Howard  et  al.  [2007, 
section  3.1]: 

1/R  =  sin  a  cots  4-  cos  cv  (1) 

where  R  is  the  distance  of  the  measured  point  from  the 
Sun  in  AU,  e  is  the  elongation  angle,  and  a  is  the  angle 
subtended  by  the  measured  point  at  the  Sun.  In  terms 
of  the  heliocentric-ecliptic  colatitude  6  and  longitude  4>, 

cosa  =  sin  0cos<t>.  (2) 

The  6  and  4>  coordinates  were  obtained  from  the 
heliographic  location  of  the  surface  feature  that  CDAW 
WG1  associated  with  each  LASCO  CME,  flare,  active 
region,  etc.,  as  given  in  the  Source  Location  column  of 
Table  2. 

[67]  On  each  of  the  (3-D)  distance-time  plots  for  the 
CDAW  events  we  then  fit  the  SMEI  ICME  data  points  writh 
a  best  fit  linear  solution  and  extrapolated  that  line  to  1  AU. 
This  1  AU  crossing  time,  given  in  the  Predicted  3-D 
column  of  Table  2,  then  represents  the  "arrival  time"  of 
the  ICME  front  at  1  AU.  This  time  can  then  be  compared 
with  the  actual  shock  AT  at  1  AU  from  the  Shock  column 
of  Table  2  and  the  algebraic  differences  tabulated  as 
shown  by  Howard  et  al.  [2006,  Table  3].  In  our  case  we  also 
compared  the  SMEI  ICME  arrival  time  in  the  Predicted 
3-D  column  with  the  storm  onset  time  in  the  Dst  Date 
column.  These  times  and  their  differences  are  given  in 
Table  3.  (This  fitting  procedure  can  also  be  used  to 
estimate  a  projected  ICME  onset  time  at  0  AU,  i.e.,  at 
the  Sun,  and  the  ICME  speed.  However,  we  did  not  use 
those  values  in  this  study.) 

[6H[  The  time  differences  are  given  in  the  AT  1CME- 
Shock  and  AT  Storm-ICME  columns  of  Table  3  for  the 
shock-ICME  AT  and  storm  onset-ICME  AT,  respectively.  In 
the  AT  lCME-Shock  column,  the  plus  sign  means  that  the 
ICME  arrived  after  the  shock,  and  in  the  AT  Storm-ICME 
column  that  the  storm  onset  followed  ICME  arrival.  This 
ordering  was  adopted  on  the  basis  of  the  expected  standoff 
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Table  3.  SMEI  Data-Based  Predictions  of  ICME  Arrival  Times  at  1  AU 


CDAW 

Storm 

Shock  Arrival 
Date  and  Time 
(UT) 

Storm  Onset 
Date  and  Time 
(UT) 

SMEI  ICME 

Arrival  Date  and  Time 
(UT) 

AT  lCME-Shock 
(hour) 

AT  Storm-ICME 
(hour) 

63 

29  Mav,  1155; 

29  May,  2200 

2003 

29  May,  2134,  19  May,  2238 

+3  05,  +4.1 

+0.4,  0.4 

66 

29  May  1831 

17  Aug.,  1345 

17  Aug,  2300 

18  Aug,  0705 

+17.3 

-15.9 

67 

29  Oct.,  0558 

29  Oct,  1500 

29  Oct,  1214;  29  Oct,  -2000 

+6.25,  +14.0 

+2.75,  -5.0 

69 

20  Nov.,  0835 

20  Nov,  1100 

19  Nov,  1939;  20  Nov,  1836 

12.9,  +10.0 

+15.35,  7.6 

70 

22  Jan.,  0105 

22  Jan,  1200 

2004 

22  Jan,  -0100;  22  Jan,  1706 

0.0,  +16.0 

+11.0,  5.1 

72 

3  Apr,  0900 

3  Apr,  1800 

2  Apr,  1634;  3  Apr,  0910; 

-16.4,  0.0,  +14.65,  +19.75 

+25.4,  +8.8,  -5.65,  10.75 

73 

22  Jul,  0945 

22  Jul,  2300 

3  Apr,  2339;  4  Apr,  0445 

22  Jul,  1732;  22  Jul,  2251 

+7.75,  +13  4 

+5.5,  +0.15 

77 

7  Nov,  0156; 

7  Nov,  2200 

7  Nov,  0324 

+1.5,  -14.6 

+18.6 

78 

7  Nov.,  1759 

9  Nov.,  0925 

10  Nov,  1000 

9  Nov,  0537;  9  Nov,  1643 

-  3.8,  +7.3 

+28.4,  +17.3 

79 

17  Jan,  0715 

17  Jan,  —1800 

2005 

17  Jan,  2152;  18  Jan,  0240 

+14,6,  +19.4 

3.9,  -8.7 

80 

21  Jan,  1652 

21  Jan,  2100 

21  Jan,  0420;  22  Jan,  0427 

12.5,  +11.6 

+16.7,  -7.45 

86 

24  Aug,  0545 

24  Aug,  1000 

24  Aug,  1919 

+13.6 

9.3 

88 

11  Sep,  0100 

11  Sep,  0300 

11  Sep,  0557 

+4.95 

2.95 

14  Dec,  1357 

15  Dec,  0000 

2006 

14  Dec,  1629 

+2.55 

+7.5 

Mean 

RMS 

+5.44  (26) 

±10.48 

+3.01  (25) 

±11  95 

distance  between  an  IP  shock  and  its  ejecta  and  the  shock - 
SMEI  auroral  onset  -  storm  onset  sequence  found  above. 
Note  that  there  can  be  more  than  one  entry  per  CDAW 
storm  on  the  basis  of  multiple  shocks  (77)  or  the  tracking  of 
multiple  ICMEs  by  SMEI  prior  to  the  storm.  For  space 
weather  prediction  purposes  we  chose  not  to  interpret 
which  ICME  was  best  associated  with  the  storm,  but  to 
group  them  together  to  find  an  uncertainty  in  their  pre¬ 
dicted  ATs.  In  a  sense  this  is  a  worst  case  scenario  but  one 
that  a  forecaster  might  typically  confront  given  only  hours 
to  make  a  storm  forecast.  We  calculated  the  algebraic  mean 
of  these  time  difference  distributions  and  found  them  to  be 
+5.44  and  +3.01  hours,  respectively,  for  ICME  AT-shock  and 
storm  onset-lCME  AT.  Thus,  on  average  for  each  event,  the 
AT  of  the  shock  was  followed  5.4  hours  later  by  the  AT  of  the 
ICME,  which  was  followed  3.0  hours  later  by  storm  onset. 

[6s]  We  can  now  compare  our  results  with  those  from 
Howard  et  al.'s  [2006]  study,  which  was  based  on  the  time 
differences  (Figure  8)  between  predicted  (LASCO  CMEs 
and  SMEI  ICMEs)  and  actual  (ACE  shocks)  arrival  times  for 
15  events  associated  with  storms  having  a  range  of  peak  Dst 
values.  The  algebraic  mean  for  ICME  AT  shock  we  calcu¬ 
lated  for  that  study  was  +4.0  hours.  This  value  is  reasonably 
consistent  with  the  +5.4  hours  we  found  for  this  study. 

[7o|  A  much  debated  subject  in  space  weather  research 
is  the  accuracy  with  which  the  geoeffective  ICME  AT  can  be 
predicted  by  a  given  procedure  or  model.  This  is  usually 
given  statistically  as  a  range  or  uncertainty  in  hours  for  a 
given  set  of  ATs.  Quantitatively,  this  can  be  stated  in  terms 
of  the  standard  deviation,  or  RMS  of  a  set  of  predicted 
versus  actual  ATs.  From  Table  3,  the  RMS  values  for  the 


ICME-shock  and  storm-ICME  time  differences  were  +10.48 
and  +11.95  hours,  respectively.  (Again,  these  "uncertainties" 
are  relative  to  the  average  ICME  transit  time  of  54  hours.) 
The  ICME-shock  RMS  time  difference  from  Howard  et  al.'s 
[2006]  study  was  ±11.4  hours;  thus,  our  results  for  the  CDAW 
storms  show  an  improved  accuracy  by  about  one  hour 
versus  Howard  et  al.'s  [2006]  study.  Although  Howard  et  al. 
[2006]  used  a  different  range  of  storms,  perhaps  the  main 
reason  for  this  difference  is  that  they  converted  SMEI  ICME 
elongation  to  distance  using  the  point  P  approximation 
whereas  we  used  the  "fixed  4>"  method. 

[?i ]  Our  uncertainty  result  using  SMEI  data  to  predict 
ICME  AT,  and  therefore  storm  onset,  of  ±10.5  hours  rep¬ 
resents  an  improvement  over  similar  values  published  in 
the  literature.  These  include  (1)  — ±16  hours  by  Gopalswamy 
et  al.  [2001]  (Note  that  they  quoted  a  mean  error  of 
10.7  hours;  this  is  not  the  RMS,  which  we  estimate  as  about 
±16  hours);  (2)  ±12.0  (HAFv2),  12.6  (ICME-ICME),  and 
11.4  hours  (1CME-IP)  by  Cho  et  al.  [2003];  (3)  ±12.2  (STOA), 
11.2  (ISPM),  and  11.6  hours  (HAFv2)  by  Fry  et  al.  [2003]; 
(4)  ±14  hours  by  Schwerin  et  al.  [2005];  and  (5)  ±11.7  (STOA), 
11.0  (ISPM),  and  11.5  hours  (HAFv2)  bv  McKetma-Laivlor 
ctal.  ]2006]. 

[72\  Thus,  these  results  demonstrate  that  the  SMEI 
observations  can  be  used  to  improve  forecasting  skill  in 
predicting  the  arrival  of  potentially  geoeffective  solar  wind 
disturbances  at  Earth.  However,  our  analysis  is  limited  in 
that  it  that  it  used  only  events  known  to  be  associated  with 
intense  geomagnetic  storms  and,  therefore,  by  definition  it 
cannot  include  'correct  null'  predictions.  The  sample  size 
is  also  small.  To  properly  evaluate  forecast  skill  using 
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standard,  e.g.,  meteorological  metrics,  it  would  be  neces¬ 
sary  to  conduct  a  double-blind  test  on  a  statistically  signif¬ 
icant  sample  of  events  without  a  priori  knowledge  of  an 
associated  storm.  This  is  beyond  the  scope  of  this  study. 

4.2.2.  Attempted  Near-Real-Time  Forecasting 
With  SMEI 

[73]  SMEI  has  demonstrated  its  ability  to  track  CMEs 
from  near  the  Sun  to  Earth,  thus  providing  a  new  capability 
for  forecasting  the  occurrence  of  storms.  SMEI  was  able  to 
detect  a  number  of  geoeffective  ICMEs  at  elongations  of 
20-30°  (from  Sun  center),  or  estimated  distances  as  far 
away  as  2/3  of  the  distance  between  the  Sun  and  Earth. 
Depending  on  the  speed  of  the  ICME  front,  these  distances 
would  correspond  to  lead  times  of  10  hours  to  1-2  days. 
SMEI  was  not  designed  to  be  an  operational  mission,  so  we 
were  usually  unable  to  use  the  SMEI  data  to  attempt  near- 
real-time  forecasts.  Therefore,  all  of  the  studies  discussed  in 
this  paper  were  retrospective  in  nature. 

[74]  However,  because  of  improved  circumstances  such 
as  an  improved  data  latency  later  in  the  mission,  we  were 
able  to  produce  near-real-time  predictions  using  LASCO 
and  SMEI  data  for  two  storms.  The  first  was  CDAW  77  in 
early  November  2004  (see  Figures  5  and  6).  In  an  AFRL 
email  issued  on  7  November  at  2300  UT,  one  of  us,  DRM, 
predicted  the  Earth  arrival  of  the  ICME  observed  by  SMEI 
on  6-7  November  to  be  '•'42:00  on  7  November.  Quicklook 
SMEI  images  also  showed  brightening  aurora  at  1400  UT, 
indicating  to  DRM  imminent  storm  onset  even  before  the 
written  email  alert.  The  ICME  was  correctly  attributed  to 
the  first  of  the  two  LASCO  halo  CMEs  observed  on  4  No¬ 
vember  (Table  2).  In  an  email  two  hours  later  we  noted  that 
the  ACE  data  indicated  "a  possible  shock  and  ejecta  and  a 
storm  with  Kp  =  7  late  on  7  November''  with  storm  onset  at 
2200  UT.  This  was  about  the  time  of  the  email  prediction; 
when  first  observed  the  ICME  was  already  at  —80°  e,  pro¬ 
viding  little  advance  warning  time.  However,  we  note 
thatour  predicted  AT  on  thebasis  of  the  current  d- 1  plot 
(Figure  6)  was  even  earlier,  0324  UT,  on  7  November. 

[75]  SMEI  observed  and  tracked  two  eastern  ICMEs  on  2 
and  4  January  2008  associated  with  a  solar  east  limb  active 
region  and  CMEs  observed  by  LASCO  on  31  December 
and  2  January,  respectively.  One  of  us,  TAK,  used  the 
SMEI  t-time  plots  to  predict  their  1  AU  crossing  times 
to  be  3  January  at  1500  UT  and  5  January  at  0030  UT,  in 
emails  on  3  January  at  1441  UT  and  4  January  at 
1536  UT.  Because  of  the  ICMEs'  east  limb  sources,  we 
did  not  expect  storms  to  occur  at  Earth  and  indeed  there 
were  none.  Thus,  these  constitute  "null"  predictions. 
There  was  possible  ICME  material  observed  in  the  ACE 
data  embedded  in  a  sector  boundary  crossing  early  on 
5  January,  in  agreement  with  the  SMEI  predicted  AT. 

5.  Discussion:  A  New  Capability  for  Forecasting 
Space  Weather 

|7f>I  SMEI  began  routine  observations  in  February7  2003, 
and  has  maintained  an  average  dutv  cycle  of  85" o  to  the 


present.  This  period  covers  the  CDAW  storms  63  -  88  occur¬ 
ring  from  May  2003  to  October  2005.  There  was  a  long  outage 
of  SMEI  data  from  28  April  to  8  June  2005,  resulting  in  a  lack 
of  SMEI  data  for  CDAW  storms  81-84.  However,  SMEI  was 
able  to  observe  the  associated  ICMEs  for  most  of  the  other 
CDAW  intense  storms.  Wre  described  the  SMEI  observations 
and  analyses  for  the  26  CDAVV  storms,  numbers  63-88,  that 
occurred  during  the  SMEI  observations  from  May  2003  to 
October  2005  as  well  as  for  the  two  other  intense  storms  that 
occurred  after  that  period  through  the  end  of  2007,  for  a  total 
of  28  storms  considered. 

[77]  The  main  statistical  results  of  this  study  are  sum¬ 
marized  in  Appendix  A.  Wre  studied  the  18  CDAW  storms 
with  adequate  SMEI  coverage  and  for  which  CDAW  WG1 
identified  the  solar  and  IP  sources.  Of  these  18  storms,  12 
were  preceded  by  a  definitely  associated  SMEI  ICME,  four 
by  a  likely  associated  ICME,  and  two  had  no  associated 
ICME.  Thus,  we  conclude  that  SMEI  likely  detected  and 
tracked  89%  (16/18)  of  the  ICMEs  that  caused  these  storms. 

[78]  The  typical  pattern  for  these  intense  events  was  IP 
shock,  SMEI  aurora  onset,  storm  onset,  with  each  interval 
about  four  hours  in  length.  The  mean  time  difference 
between  when  SMEI  first  observed  the  ICME  and  the  SMEI 
aurora  onset  time  was  23.2  hours  (range  7-36  hours).  The 
mean  time  difference  between  when  SMEI  first  observed 
the  ICME  and  the  storm  onset  time  was  27.5  hours  (range 
14-44  hours).  Therefore,  on  average  SMEI  is  capable  of  the 
first  detection  of  a  geoeffective  ICME  — 1  day  in  advance, 
yielding  a  prediction  lead  time  of —18  hours  but  with  a  large 
range.  Since  a  forecaster  requires  about  6  hours  of  SMEI  data 
to  input  into  an  arrival  time  prediction  program,  the  resul¬ 
tantwarning  time  is  reduced  to  an  average  of  —18  hours.  Of 
course,  the  time  differences  ranged  over  a  factor  of  five,  so 
the  actual  lead  time  will  be  better  or  worse  than  this, 
depending  on  the  speed  and  brightness  of  each  event. 

[79]  Using  our  method  for  converting  elongation  to  dis¬ 
tance,  we  were  able  to  make  a  prediction  of  the  ICME 
arrival  time  at  1  AU  for  each  event  and  to  calculate  the  time 
differences  for  tw  o  sets  of  data:  the  time  between  shock  and 
ICME  arrival  times,  and  the  time  between  the  ICME  arrival 
time  and  the  storm  onset  On  average  for  each  event,  the  AT  of 
the  shock  was  followed  5.4  hours  later  by  the  AT  of  the  ICME, 
which  was  followed  3.0  hours  later  by  storm  onset.  The 
algebraic  mean  for  the  ICME  AT  shock  time  differences  from 
Howard  et  al.'s  [2006J  study  was  +4.0  hours,  which  was  rea¬ 
sonably  consistent  with  the  +5.4  hours  we  found  for  this  study. 

[so]  These  results  suggest  the  following  basic  scheme  for 
using  SMEI-type  observations  to  make  storm  forecasts. 
When  a  solar  eruption  occurs  on  the  frontside,  especially  if 
within  ~45  heliodegrees  of  Sun  center,  a  forecaster  should 
check  if  an  associated  CME  is  observed  by  a  coronagraph 
within  a  reasonable  time  window.  If  so,  and  especially  if 
the  CME  is  a  partial  or  full  halo,  the  forecaster  can  then 
check  the  SMEI  observations  a  day  or  so  later,  depending 
on  the  apparent  speed  of  the  initial  CME,  for  any  possibly 
associated  ICMEs.  If  one  is  found,  then  the  leading  edge  is 
measured  over  several  successive  frames  and  an  elongation 
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versus  time  plot  constructed.  For  SMEI  we  recommend  that 
a  minimum  of  three  data  points,  thus  over  three  orbits  or 
—6  hours,  are  necessary  for  such  a  plot.  Then,  assuming  the 
solar  source  location  and  onset  time  of  the  CME  and  that  it 
propagates  radially,  a  distance-time  relation  can  be  calcu¬ 
lated  and  extrapolated  to  the  1  AU  distance.  The  storm 
onset  time  will  then  be  predicted  to  be  —3  hours  after  this 
1  AU  arrival  time.  The  prediction  program  can  be  rerun  and 
the  predicted  1  AU  AT  updated  as  more  e-t  data  points 
become  available  to  the  forecaster.  Note  that  actual  e-t  and 
d-t  plots  are  not  needed,  but  they  will  be  helpful  in  con¬ 
firming  the  predicted  AT  with  solar  wind  and  geomagnetic 
data  such  as  we  have  shown  in  this  paper. 

[si]  Of  course,  there  are  some  limitations  to  making 
such  forecasts.  If  no  associated  ICME  is  detected  by  the 
imager,  then  no  useful  forecast  can  be  made.  This  could 
lead  to  a  missed  forecast  if  there  is  a  geoecffective  ICME 
that  is  not  detected  (or  not  detected  in  time)  by  the 
imager.  In  addition,  a  SMEI-type  imager  cannot  measure 
magnetic  fields  and,  therefore,  even  if  an  ICME  is 
detected,  the  strength  and  orientation  of  the  magnetic 
field  within  or  ahead  of  it  would  remain  unknown.  How¬ 
ever,  this  is  more  likely  to  cause  uncertainty  in  the  intensity 
of  a  storm,  rather  than  no  storm  at  all,  because  most  ICMEs 
tend  to  have  geoeffective  southward  fields  somewhere 
ahead  of  or  within  themselves.  We  also  note  that  our 
analysis  was  limited  in  that  the  sample  size  was  small  and 
its  analysis  was  in  only  one  direction,  from  intense  storms 
back  to  their  solar  sources. 

[82]  The  RMS  range  between  model  or  procedure  pre¬ 
dicted  and  actual  geoeffective  ICME  ATs  is  often  used  as  a 
measure  of  the  uncertainty  in  the  procedure  or  model.  We 
calculated  the  RMS  values  for  the  ICME-shock  and  storm- 
ICME  time  differences  of  our  study  which  w  ere  ±10.48  and 
±11.95  hours,  respectively  (relative  to  the  average  ICME 
transit  time  of  54  hours).  The  ICME-shock  RMS  time  dif¬ 
ference  from  Howard  et  al.'s  [2006]  study  was  ±11.4  hours; 
thus,  our  results  for  the  CDAW  storms  show  an  improved 
accuracy  by  about  1  hour  compared  to  Howard  et  al.'s 
[2006]  study  and  other  similar  studies  in  the  literature. 
Although  one  can  argue  that  one  hour  is  of  the  order  of 
the  observational  cadence  of  both  ACE  and  SMEI,  our  re¬ 
sults  remain  comparable  to  other  studies  and,  therefore, 
demonstrate  an  improvement.  We  believe  that  much  of 
this  improvement  is  because  we  had  available  to  us  the 
CDAW  source  region  information  for  each  event  that 
allowed  us  to  convert  the  SMEI  ICME  elongations  to  dis¬ 
tances  using  the  ''fixed  4>"  method.  Howard  et  al  [2006]  used 
the  "point  P"  method.  We  note  that  with  SMEI  imagery 
alone  we  do  not  know  which  part(s)  of  the  shock/sheath/ 
ejecta  ensemble  we  are  observing  in  any  given  event  or 
whether  solar  wind  material  accumulated  with  a  snowplow- 
type  of  action  plays  a  role.  In  addition,  the  methods  in  the 
other  studies  differed  somewhat  in  whether  only  shocks,  or 
shocks  versus  ICMEs,  or  CMEs  versus  ICMEs  were  com¬ 
pared  for  arrival  times.  However,  given  the  large  distances 
and  travel  times  involved,  we  do  not  think  that  the  few 


hours  differences  between  shock  and  ejecta  arrival  times 
significantly  biases  our  result  compared  to  the  others. 

[83]  Interplanetary  coronal  mass  ejections  (ICMEs)  are  a 
primary  cause  of  severe  space  weather  at  Earth  because 
they  drive  shocks  and  trigger  geomagnetic  storms  that  can 
damage  spacecraft  and  ground-based  systems.  The  SMEI 
experiment  has  detected  geoeffective  ICMEs  at  —1/3  of  the 
distance  from  the  Sun  to  Earth  and  beyond,  corresponding 
to  an  average  advance  warning  time  of  — 1  day  (relative  to 
the  average  ICME  transit  time  of  2.2  days).  SMEI's  images 
of  ICMEs,  obtained  with  a  cadence  of  —100  min,  provide  a 
dynamic  view  of  ICME  morphology  along  trajectories 
aimed  toward  Earth.  Our  investigations  of  ICME  intensity, 
structure  evolution  and  kinematics  should  improve  our 
ability  to  forecast  storm  effects  at  Earth. 

[84]  Detecting  and  tracking  ICMEs  in  this  distance  range 
is  a  new  capability  that,  when  combined  with  other  space 
environment  sensors  and  physics-based  heliospheric 
models,  promises  to  greatly  enhance  our  ability  to  forecast 
and  study  space  weather.  The  SMEI  allows  near-real-time 
tracking  of  ICMEs  as  they  propagate  outward  from  the  Sun 
through  the  heliosphere.  Above  we  outlined  a  very  simple 
approach,  using  only  solar  and  SMEI  data,  which  a  fore¬ 
caster  might  use  to  predict  ICME  arrival  at  1  AU  and  storm 
onset  at  Earth.  STEREO  is  now  providing  a  similar  capa¬ 
bility  from  two  unique  vantage  points. 

Appendix  A:  Statistical  Results  of  SMEI-CDAW 
Storm  Study 

[85]  The  main  statistical  results  from  study  of  the  CDAW 
storms  that  involved  SMEI  data  are  listed  here  and  sum¬ 
marized  in  section  5. 

Al.  CDAW  Storm:  ICME  Associations 

[86]  Of  the  28  CDAW  storms  studied,  7  had  no  or 
inadequate  SMEI  data,  and  3  more  were  identified  as 
CIR  driven.  Of  the  remaining  18  storms  with  adequate 
SMEI  coverage  and  for  which  CDAWT  WG1  identified  the 
solar  and  IP  sources,  12  were  preceded  by  a  definitely 
associated  SMEI  ICME,  four  by  a  likely  associated  ICME, 
and  two  had  no  associated  ICME.  In  addition,  for  these  18 
storms,  9  of  the  LASCO  CMEs  and  SMEI  ICMEs  were 
likely,  5  possibly,  3  questionably,  and  1  unlikely  associated. 

A2.  SMEI  ICME  Spans  and  Elongations 

[87]  The  ICME  spans  were  measured  from  the  SMEI  data 
and  are  presented  in  Table  2.  The  mean  value  was  78.63°  and 
the  range  was  30-13°.  There  were  11  measured  values  and  for 
an  additional  5  we  could  only  measure  lower  limits.  The 
elongation  angles  of  the  ICMEs  at  the  time  of  the  onset  of  each 
storm  are  presented  in  Table  2.  The  mean  value  was  78.67° 
and  the  range  was  45-105°.  There  were  15  measured  values. 

A3.  Timing  Results 

1 88 1  The  observed  time  differences,  .AT,  between  the 
onset  time  of  each  storm  and  the  onset  time  of  the  SMEI 
aurora  are  listed  in  Table  2  (column  f^lff  St-Au).  The  mean 
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value  was  +4.085  hours,  and  the  range  was  -2.5-8.15 
hours.  There  were  10  measured  values.  The  observed  time 
differences,  AT,  between  the  onset  time  of  each  SMEI 
aurora  and  the  shock  arrival  time  at  1  AU  are  listed  in 
Table  2,  column  19.  The  mean  value  was  +3.88  hours,  and 
the  range  was  -22.9  to  +15.4  hours.  There  were  14  mea¬ 
sured  values  and  for  an  additional  4  we  could  only 
measure  lower  limits.  There  were  only  three  events  in 
which  the  shock  preceded  the  aurora  onset. 

[89]  The  typical  sequence  of  these  three  key  events  was 
the  arrival  at  1  AU  of  the  IP  shock,  onset  of  the  SMEI 
aurora,  and  onset  of  the  storm.  From  above,  each  of  these 
time  intervals  was  about  four  hours  in  duration.  The  mean 
time  difference,  AT,  between  when  SMEI  first  observed  the 
ICME  and  the  onset  of  the  SMEI  aurora  was  23.23  hours, 
with  a  range  of  7-36  hours  (Table  2  column  tdiff  Au-FO). 
There  were  11  measured  values.  The  mean  time  difference 
between  when  SMEI  first  observed  the  ICME  and  the  storm 
onset  time  was  27.53  hours,  with  a  range  of  14-44  hours 
(Table  2  column  Au-FO).  There  were  15  measured  values. 

[90]  Using  our  method  for  converting  measured  SMEI 
elongation  angles  to  distance,  we  could  predict  the  ICME 
arrival  time  at  1  AU  for  each  event  and  use  it  to  calculate 
the  time  differences  for  two  sets  of  data:  the  time  between 
shock  and  ICME  arrival  times  (ATs),  and  the  time  between 
the  ICME  arrival  time  and  the  storm  onset.  These  results 
are  given  in  Table  3.  The  mean  value  between  the  AT  of 
the  shock  and  the  AT  of  the  ICME  was  +5.44  hours.  The 
mean  value  between  the  storm  onset  and  the  AT  of  the 
ICME  was  +3.01  hours.  Thus,  the  typical  sequence  began 
with  the  arrival  of  the  shock,  followed  5.4  hours  later  by 
the  predicted  arrival  of  the  ICME,  followed  3.0  hours  later 
by  the  onset  of  the  storm. 

[91 ]  The  RMS  range  between  model-  or  procedure- 
predicted  and  actual  geoeffective  ICME  ATs  is  often  used 
as  a  measure  of  the  uncertainty  in  the  procedure  or  model. 
The  RMS  value  for  the  ICME-shock  time  differences  was 
±10.48  hours,  with  26  values.  The  RMS  value  for  the  storm- 
ICME  time  differences  was  ±11.95  hours,  with  25  values. 
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